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PROJECT  OVERVIEW 


OBJECTIVES 

The  general  objective  of  the  project  is  fundamental  mathematical  modeling  of  a  complex  TSA 
system  with  electrothermal  desorption  step,  with  adsorbers  assembled  of  one  or  more 
activated  carbon  fiber  clot  (ACFC)  cartridge-type,  radial  flow  fixed  beds  and  with  in- vessel 
condensation  of  the  desorbed  vapor.  During  the  first  stage  of  the  project,  a  mathematical 
model  of  a  single  cartridge  without  condensation  was  developed  and  used  for  simulation. 

The  objective  of  this,  second  stage  of  the  project,  is  to  model  the  complete  TSA  system,  for 
adsorbers  with  one  and  more  cartridges.  Besides  the  material  and  energy  balances,  these 
models  have  to  include  the  momentum  balances,  in  order  to  define  the  complex  gas  flow  in 
the  complex  geometry  of  the  adsorber.  The  models  have  to  describe  adsorption, 
electrothermal  desorption  and  electrothermal  desorption  accompanied  with  condensation  of 
the  desorbed  vapor,  as  well  as  the  complete  TSA  cycle.  These  models  will  be  used  for 
prediction  of  the  velocity,  pressure,  concentration  and  temperature  profiles  in  the  adsorbers 
and  calculation  of  the  mass  of  condensed  liquid  and  used  electric  energy.  Applicability  of  the 
developed  mathematical  models  for  optimization  of  the  TSA  process,  with  regard  to  the 
separation  and  purification  factor,  regeneration  of  the  condensed  vapor  and  used  energy,  will 
also  be  investigated. 


STATUS  OF  EFFORT 

Models  for  two  types  of  adsorbers,  one  with  only  one,  and  the  other  with  two  cartridges,  have 
been  developed  in  this  stage.  These  adsorber  types  will  be  called  1 -cartridge  and  2-cartridges 
adsorber,  respectively,  through  this  manuscript.  For  each  adsorber  type,  three  models  were 
built,  in  order  to  describe  three  stages  of  a  compete  TSA  cycle:  adsorption,  electrothermal 
desorption  and  electrothermal  desorption  with  in-vessel  condensation.  These  models  were 
built  using  Femlab,  a  specialized  software  tool  for  modeling  of  complex  systems  with 
complex  geometry.  In  order  to  describe  the  complete  TSA  cycle,  the  models  for  the  three 
stages  were  integrated,  by  using  a  combination  of  Femlab  and  Matlab.  The  models  were 
successfully  used  for  simulation  of  separate  stages  of  the  process  and  of  the  complete  TSA 
cycles,  as  well  as  for  their  optimization. 


ACCOMPLISHMENTS  /NEW  FINDINGS 

The  accomplishments  of  this  project  can  be  perceived  from  two  aspects: 

1 .  The  developed  models  represent  a  mathematical  tool  which  can  be  used  for  simulation 
of  the  investigated  TSA  system,  its  analysis,  predicting  of  the  concentration  and 
temperature  profiles  in  the  system  and  its  optimization.  This  was  the  main  objective  of 
this  project. 

2.  By  developing  Femlab  and  combined  Femlab/Matlab  models  of  this,  very  complex 
system,  involving  a  number  of  coupled  processes  and  phenomena  taking  place  in  a 
complex  geometry,  a  progress  in  the  Femlab  and  Femlab/Matlab  modeling  has  been 
achieved. 
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DETAILED  REPORT 


1.  Introduction 

The  idea  about  regeneration  of  used  adsorbents  by  direct  heating  of  the  adsorbent  particles  by  passing 
electric  current  through  them  (Joule  effect),  was  first  published  in  the  1970s  [1].  Desorption  process 
based  on  this  principle  was  later  named  electrothermal  desorption  [2].  It  was  recognized  as  a 
prospective  way  to  perform  desorption  steps  of  TSA  cycles.  At  the  same  time,  fibrous  activated  carbon 
was  recognized  as  a  very  convenient  adsorbent  form  for  its  realization.  Electrothermal  desorption  has 
some  advantages  over  conventional  methods,  regarding  adsorption  kinetics  and  dynamics  [3,  4]  and 
energy  efficiency  [2,  5].  Some  industrial  applications  of  electrothermal  desorption  have  been  reported 
recently  [6,  7].  Nevertheless,  the  development  of  mathematical  description  of  processes  based  on 
electrothermal  desorption  hasn’t  been  following  the  development  of  the  process  itself,  so  far. 

A  new  TSA  process  with  electrothermal  desorption  step,  based  on  adsorbers  assembled  of  one  or  more 
annular,  cartridge-type,  fixed-beds,  with  in- vessel  condensation,  has  been  presented  recently  [5,  8,  9]. 
Its  schematic  representation  is  given  in  Figure  1. 


EfFlu  ent 


Condensate 


Figure  1.  Overall  schematic  representation  of  the  ACFC  adsorption  -  rapid  electrothermal  desorption  system  (from  Ref.  [5]) 


The  aim  of  this  project  is  to  develop  a  rigorous,  fundamental  mathematical  model  of  this  system, 
which  could  be  used  for  simulation,  analysis  and  optimization  of  the  TSA  process  developed  by 
Sullivan  [5]. 


1.1.  Description  of  the  system 

The  adsorbers  of  the  TSA  system  shown  in  Figure  1  are 
composed  of  one  or  more  annular,  radial-flow  cartridge-type 
adsorbent  beds.  A  single  cartridge  is  schematically  shown  in 
Figure  2.  It  is  formed  as  a  cylindrical  roll  of  activated  carbon 
fiber  cloth  (ACFC),  spirally  coiled  around  a  porous  central 
pipe.  The  gas  flow  through  the  adsorber  is  in  the  radial 
direction.  During  the  desorption  step,  electric  current  is 
passed  through  the  activated  carbon  cloth  in  the  axial 
direction,  causing  heat  generation,  heating  of  the  adsorbent 
and  desorption. 


GAS  OUT 


|  CAS  IN 

Figure  2.  Schematic  representation  of  the 
annular  -  radial  flow  -  ACFC  adsorbent  bed 
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The  cartridges  are  situated  in  a  cylindrical  adsorber  vessel.  During  desorption,  the  hot  gas  steam  rich 
with  the  desorbed  vapor  comes  into  contact  with  the  cold  vessel  wall  and  the  vapor  condenses.  In  that 
way,  the  adsorber  vessel  serves  as  a  passive  condenser.  The  bottom  outlet  of  the  adsorber  vessel  is 
funnel-shaped,  in  order  to  facilitate  draining  of  the  condensed  liquid. 


In  this  phase  of  the  project,  adsorbers  with  one  and  two  cartridges 
will  be  modeled.  The  1 -cartridge  adsorber  is  shown  in  Figure  3  and 
the  2-cartridges  adsorber  in  Figure  4. 


In  Figures  3  and  4  the  directions  of  the  fluid 
flow  are  also  shown,  in  order  to  clarify  the 
processes: 

-  For  the  1 -cartridge  adsorber  the  flow 
directions  during  adsorption  and  desorption 
are  the  same:  the  gas  enters  the  adsorber 
through  the  central  tube,  from  the  bottom,  and, 
after  passing  through  the  adsorbent  bed  in  the 
radial  direction,  exits  from  the  annular  space 
around  the  adsorbent  bed  through  the  funnel- 
shaped  bottom  outlet.  During  desorption,  the 
desorbed  vapor  is  partially  condensed  at  the 
cold  outer  wall  of  the  adsorber,  and  the  outlet 
stream  contains  both  gas  and  the  condensed 
liquid. 


(A)  (B) 


Gas  out  Gas  in 


t  1 


1 

Gas  and 


Figure  3.  1 -cartridge 
adsorber 


Figure  4.  2-cartridge  adsorber: 
(A)  adsorption,  (B)  desorption 


-  For  the  2-cartridges  adsorber,  the 
flow  directions  during  adsorption  and 
desorption  are  somewhat  different.  During 
adsorption  the  feed  stream  enters  the  adsorber  from  the  bottom,  through  the 
central  tube  of  the  bottom  cartridge,  and  the  purified  stream  exits  the  adsorber  at  the  top,  through  the 
central  tube  of  the  top  cartridge  (Fig.  4A).  Between  the  inlet  and  the  outlet,  the  stream  flows  through 
both  adsorbent  beds  in  the  radial  direction,  but  in  the  opposite  directions  (through  the  bottom  cartridge 
from  the  center  to  the  periphery,  and  through  the  top  cartridge  from  the  periphery  to  the  center). 
During  desorption,  two  inert  streams  are  introduced  into  the  adsorber,  from  the  bottom  and  the  top, 
through  the  central  tubes  of  both  cartridges,  and  the  gas  and  the  condensed  liquid  exit  the  adsorber 
through  the  funnel  shaped  bottom  outlet  of  the  annular  tube  around  the  cartridges  (Fig  4B).  It  should 
be  noticed  that  the  gas  cannot  flow  directly  between  the  central  tubes  of  the  top  and  the  bottom 
cartridge.  Also,  the  electric  current  supply  for  the  two  cartridges  is  separate. 


The  inlet  gas  flow-rate  during  adsorption  is  generally  several  times  higher  than  during  desorption.  The 
inlet  gas  during  desorption  is  pure  inert. 


1.2.  Brief  survey  of  the  first-year  results 

During  the  first  year,  the  modeling  efforts  were  concentrated  on  modeling  of  a  single  annular, 
cartridge-type,  fixed-bed  adsorber,  without  condensation.  Two  deterministic  models  were  postulated: 
the  first  one  assuming  uniform  adsorbent  density  throughout  the  adsorbent  bed,  and  the  second  one 
taking  into  account  the  structure  of  the  bed  made  of  layers  of  ACFC.  Both  models  were  obtained  as 
complex  sets  of  coupled  nonlinear  PDEs,  ODEs  and  algebraic  equations.  Numerical  procedures  for 
solving  the  model  equations  were  established  and  used  for  simulation  of  adsorption,  electroresistive 
heating,  electrothermal  desorption  and  consecutive  adsorption-desorption,  using  Matlab  software. 
Similar  results  were  obtained  with  the  two  models.  The  results  of  this  stage  are  presented  in  the  Final 
Performance  Report  for  project  FA8655-03-1-3010  [10]. 
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1.3.  Outline  of  the  second  year  project 

The  TSA  system  defined  in  Ref.  [5]  is  very  complex,  regarding  both  its  geometry  and  the  processes 
involved.  For  that  reason,  a  decision  was  made  to  perform  the  modeling  of  the  complete  system  by 
using  some  specialized  modeling  software.  Our  choice  was  to  use  FEMLAB,  relatively  new  modeling 
software  based  on  the  finite-element  method  for  solving  partial  differential  model  equations. 

Analysis  of  the  TSA  system  with  electrothermal  desorption  and  in-vessel  condensation  shows  that  the 
complete  cycle  can  be  broken  down  into  three  consecutive  steps: 

adsorption, 

desorption  without  condensation, 
desorption  accompanied  with  condensation. 

Although  these  three  stages  can  be  described  by  the  same  set  of  model  equations,  the  boundary 
conditions  corresponding  to  the  three  stages  are  essentially  different,  so  three  different  Femlab  models 
need  to  be  built.  In  order  to  simulate  a  complete  TSA  cycle,  these  3  models  have  to  be  run  in 
consequence  with  automatically  switching  from  one  model  to  the  next,  when  predefined  conditions  are 
reached.  The  most  convenient  way  to  achieve  that  is  to  build  a  combined  Femlab/Matlab  model  for  the 
complete  cycle.  Naturally,  separate  models  have  to  be  built  for  the  1 -cartridge  and  2-cartridges 
adsorbers. 

Taking  all  this  into  consideration,  these  are  the  main  phases  of  the  second  year  project: 

1.  Building  Femlab  models  for  adsorption,  desorption  and  desorption  accompanied  with 
condensation  for  1 -cartridge  adsorber; 

2.  Building  Femlab  models  for  adsorption,  desorption  and  desorption  accompanied  with 
condensation  for  2-cartridges  adsorber; 

3.  Building  Femlab/Matlab  models  for  integral  TSA  cycles  for  1 -cartridge  and  2-cartridges 
adsorbers; 

4.  Simulation  of  adsorption,  desorption  and  desorption-condensation  in  1 -cartridge  and  2- 
cartridges  adsorbers; 

5.  Simulation  of  complete  TSA  cycles  for  1 -cartridge  and  2-cartridges  adsorbers; 

6.  Investigation  of  the  influence  of  the  main  operational  parameters  of  the  TSA  process  and 
optimization  bases  on  it. 


2.  The  FEMLAB  models 

Six  Femlab  models  were  built: 

•  Model  A  l  -  model  for  adsorption  in  the  1-cartrdge  adsorber 

•  Model  D  l  -  model  for  desorption  without  condensation  in  the  1 -cartridge  adsorber 

•  Model  DC  l  -  model  for  desorption  with  condensation  in  the  1 -cartridge  adsorber 

•  Model_A_2  -  model  for  adsorption  in  the  2-cartrdges  adsorber 

•  Model_D_2  -  model  for  desorption  without  condensation  in  the  2-cartridges  adsorber 

•  Model_DC_2  -  model  for  desorption  with  condensation  in  the  2-cartridges  adsorber. 
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Model  assumptions: 

The  following  assumptions  were  used  in  setting  up  the  models: 

•  The  adsorbent  beds  are  treated  as  homogeneous,  with  uniform  adsorbent  porosity  and  density. 

•  The  mass  and  heat  transfer  resistances  on  the  particle  scale  are  neglected. 

•  The  fluid  phase  is  treated  as  an  ideal  gas  mixture  of  the  inert  and  the  adsorbate. 

•  All  physical  parameters  and  coefficients  are  considered  as  constants. 

•  The  electric  resistivity  of  the  ACFC  adsorbent  is  temperature  dependent.  Linear  temperature 
dependence  is  assumed,  based  on  experimental  results  reported  in  Ref.  [5]. 

•  The  electric  power  during  electrothermal  desorption  is  supplied  under  constant  voltage 
conditions. 

•  The  condensation  at  the  adsorber  wall  is  dropwise.  This  assumption  is  based  on  experimental 
observations  reported  in  Ref.  [5]. 

•  The  volume  of  the  condensed  liquid  is  neglected,  i.e.,  it  is  assumed  that  the  liquid  drops  don’t 
influence  the  gas  flow. 

•  The  heat  resistance  and  heat  capacity  of  the  adsorber  wall  are  neglected,  so  that  the  wall 
temperature  is  equal  to  the  temperature  of  the  environment. 

•  Initially,  the  adsorbate  concentrations  in  both  phases  are  in  equilibrium  and  uniform 
throughout  the  adsorbent  bed.  The  temperatures  of  both  phases  are  initially  equal  and  uniform 
throughout  the  adsorbent  bed. 


The  main  differences  from  the  models  given  in  Ref  [10] 

Owing  to  the  potentials  of  the  Femlab  software  to  model  rather  complex  systems,  some  of  the 
assumptions  used  during  the  first  phase  of  the  project  (Grant  No.  FA8655-03-1-3010  [10])  have  been 
omitted  it  the  models  presented  in  this  report.  These  are  the  basic  differences  between  the  current 
Femlab  models  and  the  models  developed  during  the  first  year  of  the  project: 

-  Two-dimensional  models  are  used  instead  of  one-dimensional  (the  gradients  both  in  the  axial 
and  in  the  radial  direction  are  taken  into  account). 

-  The  momentum  transport  equations  and  the  continuity  equations  for  the  inert  gas  are 
incorporated  in  the  models  and  solved  simultaneously  with  the  mass  and  energy  balance  equations. 

-  Instead  of  assuming  perfect  mixing  in  the  inlet  and  outlet  tubes,  the  velocity  field  is 
calculated  using  the  model  equations. 

-  The  gas  pressure  and  density  in  the  system  are  not  assumed  constant,  nor  uniform.  They  are 
calculated  from  the  model  equations. 

-  Condensation  at  the  adsorber  wall  is  taken  into  account. 


2.1.  Model  building  in  Femlab 

2.1.1.  Definition  of  the  system  geometry 

When  building  a  model  in  Femlab,  it  is  necessary  to  first  define  the  space  dimension.  Axially 
symmetrical  2D  space  was  used  for  both  adsorber  types. 

The  next  step  is  to  define  the  system  geometry.  This  is  done  by  direct  plotting  in  the  Femlab  window. 
The  system  geometry  is  shown  in  Figure  5:  Figure  5a  for  the  1 -cartridge  adsorber  and  Figure  5b  for 
the  2-cartridges  adsorber.  Figure  5  represents  the  halfs  of  the  axial  cross-sections  of  the  adsorbers 
(axial  symmetry  2D  models  are  used).  The  main  dimensions  of  the  adsorbers  are  also  shown. 
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Figure  5.  Definition  of  the  model  geometry,  with  the 
main  dimensions:  (a)  for  the  1 -cartridge  and  (b)  for  the 
2-cartridge  adsorber 


2.1.2.  Definition  of  the  model  equations 

The  next  step  is  to  define  the  model  equations,  by 
choosing  the  appropriate  application  modes.  Figure  6 
shows  the  Model  Navigator  window  which  is  used  for 
defining  the  space  dimension  and  selecting  the 
application  modes. 

The  following  application  modes  were  used  in  our 
Femlab  models: 

1 .  N on-isothermal  flow  -  for  defining  the  gas  flow 
in  the  central  and  annular  tubes; 

2.  Brinkman  Equations  -  for  defining  the 
momentum  balance  for  the  packed  adsorbent 
bed(s); 

3.  Convection  and  Diffusion  -  for  defining  the 
mass  balance  for  the  gas  phase; 

4.  Diffusion  -  for  defining  the  mass  balance  for 
the  solid  phase; 

5.  Convection  and  Conduction  -  for  defining  the 
heat  balance  for  the  gas  phase; 


6.  Heat  transfer  by  Conduction  -  for  defining  the  heat  balance  for  the  solid  phase; 

7.  Conductive  Media  DC  -  for  defining  heat  generation  by  Joule  effect. 


Although  the  Coductive  Media  DC 
application  mode  is  not  needed  in  the 
adsorption  models  (there  is  no  electrical 
heating  during  adsorption),  all  7  application 
modes  were  used  in  all  models,  in  order  to 
enable  their  integration  into  the  models  of 
complete  TSA  cycles. 

The  geometry  of  the  adsorbers  assumes 
definition  of  their  parts,  i.e.  subdomains. 

Three  subdomains  are  defined  for  the  1- 
cartridge  adsorber  (see  Fig.  5a): 

Rl-  the  central  tube 
R2-  the  adsorbent  bed 
COl-  the  annular  space  around 
the  adsorbent  bed. 

On  the  other  hand,  5  subdomains  are  defined 
for  the  2-cartridges  adsorber  (Fig.  5b): 

Rl  -  the  central  tube  of  the 
bottom  cartridge 
R2  -  the  adsorbent  bed  of  the  bottom  cartridge 
R3  -  the  central  tube  of  the  top  cartridge 
R4  -  the  adsorbent  bed  of  the  top  cartridge 

COl  -  the  annular  space  between  the  cartridges  and  the  outer  adsorber  wall 


Figure  6.  The  Model  Navigator  window 


By  selecting  the  application  modes  one  by  one  in  the  Multiphysics  menu,  it  is  possible  to  access  the 
underlying  PDEs  and  define  their  coefficients  and  initial  conditions,  using  the  Subdomain  Settings 
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Menu  (an  example  is  shown  in  Figure  7),  and  the  boundary  conditions,  using  the  Boundary  Settings 
Menu  (an  example  is  shown  in  Figure  8),  for  each  subdomain. 

Checking  “Active  in  this  domain”  activates  the  chosen  application  mode  in  the  corresponding 
subdomain.  An  overview  of  the  application  modes  and  the  subdomains  in  which  they  are  active  is 
given  in  Table  1. 


Table  1.  An  overview  of  the  application  modes  and  the  subdomains  in  which  they  are  active 


Application  mode 

1-catridge  adsorber 

2-catridge  adsorber 

R1 

R2 

COl 

R1 

R2 

R3 

R4 

COl 

Non-isothermal  flow 

+ 

+ 

+ 

+ 

+ 

Brinkman  Equations 

+ 

+ 

+ 

Convection  and  Diffusion 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

Diffusion 

+ 

+ 

+ 

Convection  and  Conduction 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

Heat  transfer  by  Conduction 

+ 

+ 

+ 

Conductive  Media  DC 

+ 

+ 

+ 

Owing  to  the  complexity  of  our  model,  with  a  number  of  coupled  nonlinear  differential  equations  and 
interlinked  variables,  the  default  Femlab  PDEs  had  to  be  modified.  The  Femlab  window  Subdomain 
Settings  -  Equation  System ,  such  as  the  one  shown  in  Figure  9,  is  used  for  this  purpose. 


Equation 

Tgit  =  Tgit0 


Boundary  selection 


il 

nil 

gj 

i 

|4 

1 

is 

is 

_ 

_ U 

□  Select  by  group 
0  Interior  boundaries 


boundary  conditions 

Boundary  condition: 

Temperature 

1 

Quantity 

Value.'Expression 

Description 

% 

15 

Inward  heat  flux 

Tgit0 

293.15 

Temperature 

|  OK  1  [  Cancel  \  [  Apply  ] 


Figure  8.  Boundary  Settings  Menu 


Figure  7.  Subdomain  Settings  Menu 
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Figure  9.  Subdomain  Settings  -  Equation  System  Menu 


30 

After  modifying  the  PDEs  and  defining  the 
constants  and  the  parameters  in  the  system,  the 
mesh  for  the  finite  element  method  is 
generated.  Figure  10  shows  the  initial  meshes  2D 
for  our  models  (Fig  10a  for  the  1 -cartridge  and  ^ 
Fig  10b  for  the  2-cartridges  adsorber). 

The  models  were  solved  using  the  Direct 
(UMFPACK)  solver. 

ID 


5 


0 


■5  0  5  -10  0  10 

r  [cm]  r  fcm] 

Figure  10.  The  initial  mesh  for  the  finite-difference  method 
for:  (a)  1 -cartridge  and  (b)  2-cartridges  models 


(a)  (b) 


2.2.  The  underlying  equations 

The  following  equations  underlay  the  Femlab  models  presented  in  the  previous  section: 

1)  Momentum  balances  and  continuity  equations  for  the  inert  gas: 

For  the  central  and  outer  annular  tubes : 


8u 

nr, — + v 


g  dt 


du 

-2r|a- — 
or 


r  du  dv'' 

f 

f  / 

- 1 - 

=  — 

r 

Pc 

V  dz  dr  ) 

V 

r  g 

V  v 

du  du 

u - h  v- 

dr  dz 


)+dp" 

\ 

+  2[i- 

J  dry 

r) 

dv 

rp’T,+v 


' dv  du ' 

.  dv 

(  f 

-r\i 

Kdr  dz  y 

~2r^ 

dz 

=  —r 

Pg 

K  v 

dv  dv 

u - h  v — 

dr  dz 


+ 


dp 

dz 


(1) 

(2) 
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(3) 


f 

f 

(  du 

dv^ 

\ 

Per 

r 

— + — 

+u 

rg\ 

V 

V 

Kdr 

dz  J 

J 

■2^u+2^ 

dr  dz 


=  0 


For  the  adsorbent  bed(s): 


dii 


8  dt 


du 

-2r\i—  + 
or 


du  SvY 

- 1 - 

^  dz  dr  )  J 


A 


p  dp 

—  U  H — — 

k  dr 


A  +  2(i-A 

r 


dv 

rPg^  +  V 


8  dt 


-m 


dv  du 


^ dr  dz  J 


■2 

dz 


Li  dp 

— v  +  — 
k  dz 


f 

f 

(  du 

dv^ 

\ 

Pa 

r 

- 1 - 

|  +  u 

K  g 

V 

V 

ydr 

dz  ) 

J 

dpg  dp 

+  — -u  +  — - 
dr  dz 


=  0 


(4) 

(5) 

(6) 


The  following  notations  are  used  in  these  equations  :  u  and  v  -  gas  velocities  in  the  radial  and  in  the 
axial  direction,  respectively,  p  -  pressure,  pg  -  gas  phase  density,  t  -  time,  r  -  radial  coordinate,  z  - 
axial  coordinate,  adsorbent  bed  permeability  and  p  -  dynamic  viscosity  of  the  carrier  gas  (inert). 

The  adsorbent  bed  permeability  was  calculated  using  the  following  equation: 


k  =  -p— 

Vp 

The  gas  density  was  calculated  as: 
P 

Pg  RgTg(\  +  C ) 


(7) 


(8) 


where  Tg  and  C  are  the  temperature  and  adsorbate  concentration  in  the  gas  phase  and  Rg  the  universal 
gas  constant. 

2)  Heat  balances  for  the  gas  phase: 

For  the  central  tube(s): 


td_ 
dt 1 


£7, 

u  dr 


-rDt‘ 


2k 

li<  dz 


(  dC  dC\ 
s~pv~g{U  dr+V  dz  , 


=  -rpepX 


dF  dT 

-  rPg  (cPg  +  cpvC>  ~  rPg  ( cPg  +  cpvc>-^- 


For  the  adsorbent  bed(s): 

d 


+c-c)r«i+v  -rD?h2--rD‘i  -£ 


2k 

b  dr 


2k 

dz 


f 


(  dC  dC 


y  V  dr  dz 

For  the  outer  annular  tube : 


dT 


dT 


KatTs~T  )-p  c  T  u  —  +  v—  —  rp  (c  +  c  C)u  — - rp  (c  +c  C)v— 


dr 


dz 


_e 

dt 


k  ]+ vf-  '•o?  L  gt  -  rDi‘  |  ^ ' 


dC  dC 


s  pv  8{  dr  dz 


lot  Qr 

dT 


l°t  dz 


~rPgcpvTg  u  —  +  v—  - rpg (cpg  +  cpvC)u  — - rps(cps+cpvC)v 


dr 


g  y  pg  pv 


2k 

dz 


(9) 


(10) 


(11) 
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3)  Adsorbate  balance  for  the  gas  phase: 

For  the  central  tube(s): 


9C  Y7 

r - +  V| 

dt 


Dmr  it  dC 

^mz  it 

dE 

( 

pg  dr 

p« 

dz 

J 

V 

8C  8C ^ 
ru - h  rv 


dr 


dz 


For  the  adsorbent  bed(s)\ 


dC  V7 

r - h  V 

dt 


f 


Dmr\b  dC  A 


-  r- 


V 


Pg  dr  pg 
For  the  outer  annular  tube : 


-r- 


mz  \b 


dC 


dz 


=  r^-(C* -C)d 


dC 


g 


ru - brv  — 

V  dr  dz 


dC' 

J 


dC  V7 

r - h  V 

dt 


f 


Dmr\ot  dC  Dm 


V 


dr 


Lee' 

l  r 

r  &J 

! 

1 

dC  dC \ 

ru - h  rv - 

dr  dz  J 


(12) 


(13) 


(14) 


4)  Electric  current  balance  for  resistive  heating: 

f  1  CST  T  1  CST  t\ 


1  dU  ]_dU 

p  dr  p  dz  j 


■  0 


(15) 


In  this  equation  U  is  the  electric  potential  and  p  the  electric  resistivity  of  the  adsorbent,  which  is 
temperature  dependent.  Linear  temperature  dependence,  which  was  obtained  experimentally  for 
ACFC  adsorbent  [5],  was  used  in  our  calculations: 


p  =  p0(l+AX-A)) 


(16) 


5)  Heat  balance  for  the  solid  phase  within  the  adsorbent  bed(s): 


>Vh  f[(^  +  cplq)Ts\- v( rD';:  ^  +  rD, hs  ^  ^  -  -  8Qd 


dr 


dT 


dz 


=  r 


dq 


,v  +  rP b^ads  —  -  rhba(Ts  -  T  ) 
dV  dt 


(17) 


where  Ts  is  the  solid  temperature  and  (5 Qel  /  dV)  the  electric  power  supply  per  unite  volume  of  the 
adsorbent  bed: 


SOL  =  I 

dV  p 


dU 
\  dr  j 


+ 


dU 
ydz  j 


1 


p0(l  +  6(7;-7i)) 


f  a tt\ 2  f 


dU 

V  Sr  j 


+ 


dU 
v  dz  J 


(18) 


6)  Adsorbate  balance  for  the  solid  phase  within  the  adsorbent  bed(s): 

p  b^-  =  kma(C-C*)  (19) 

dt 

where  q  is  the  adsorbate  concentration  in  the  solid  phase  and  C*  the  adsorbate  concentration  in  the  gas 
phase  in  equilibrium  with  the  solid  phase.  The  Dubinin-Radushkevich  equation  was  used  to  describe 
the  adsorption  equilibrium: 


f 

W  =  W0ex  p 

V 


RT 


g  s  ln^— 


Pajj 


(20) 


In  this  equation  W  is  the  volume  of  the  adsorbed  phase  per  unite  mass  of  adsorbent,  W0-  the  volume  of 
the  micropores  per  unit  mass  of  adsorbent,  E-  the  adsorbate  energy  of  adsorption,  p°- the  adsorbate 
saturation  pressure,  pA  -  the  adsorbate  partial  pressure.  The  saturation  pressure  was  calculated  using 
the  Wagner  equation: 
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In 


[p°] 

'  VPAx  +  VPBx  5  +  VPCX 3  +  VPDx6  ' 

UcJ 

l  J 

(21) 


where  x  =  \-(Ts/Tc) ,  pc  and  Tc  are  the  critical  pressure  and  temperature  and  VP  a,  VPb,  VPc  and  VPD 
are  Wagner  constants. 

Based  on  equation  (20),  the  following  relation  between  the  equilibrium  adsorbate  concentration  in  the 
gas  phase  C*  and  its  concentration  in  the  solid  phase  q  and  temperature  Ts,  was  obtained: 


C 


P  exP| 


-In 


'  MaI 

KPaWo  J 


(22) 


(Ma  and  pA  are  the  adsorbate  molar  mass  and  density,  respectively). 

Boundary  conditions: 

Equations  (1-19)  are  common  for  all  three  processes  (adsorption,  electrothermal  desorption  and 
electrothermal  desorption  with  condensation)  and  for  both  adsorber  types  (1 -cartridge  and  2-cartridges 
adsorbers),  and  consequently  for  all  six  Femlab  models.  Nevertheless,  some  of  the  boundary 
conditions  for  these  6  models  are  different. 


Bl)  Common  boundary  conditions  for  all  models  for  the  1-cartridge  adsorber 

z  =  0,  r  g  (0,q) :  w  =  0,  v  =  Tg  =  Tgm,  C  =  Cin, 

r,  7ip 


r\  Wg 

z  =  H  +  h{ ,  re^/j) 

(/?!  ,hx+H):  p\it=p\b. 


:  u  =  0,  v  =  0,  -D*  1^  =  0,  -  Dmz  \  —  =  0 

hi  Qz  'it 


dz 


(23) 

(24) 


r  =  r\ ,  z 

f  QJ1  \ 

-D«\„it+v*(cn+c-’c)uT‘ 

Jit 


u\.  =  It  ,  ,  V.  =  V,, 
b 9  it  b 9 


~D-lf  *uC 


V 

^  f  ,  QC 

-DA  —  +  uC 


I  it  ™\b 

^  (  .  QT 

-D!-i-P  +  P,(crt+cp,c)uTt 


Jb 


+  hsl(l-zb)(Ts-Tg) 


nhs  A 

'D<r  AZ  =  nsAls 


Jit  v  Or  Jb 

=  hsX(Tg-Ts),  J  =  0 

0,v  =  0,-DZgld^  =  0,-Dmz\b^  =  0,U  =  U0,-D^  =  0 


(25) 


dz 


z  =  hx,  r  e  (rx,r2) : 

dT  r)( ^  r)  T 

z  -  hx+  H,  r  e  (rx,r2) :  u  =  0,  v  =  0, -D*\  =  0,  -Dmz\  —  =  0,  U  =  0,  D*  ~A  =  0 

'b  dz  D  dz  dz 

dT, 

Z  E  (0. 

+  H):  p\=p\  ,  u\  =u\  ,  v|  =v|  , 

'  b  ot’  I  b  I  ot9  I  b  \oP 

h  ,  dT 

_Dhg\  _L  +  p  rc  +c  C)uT 

tr  I °t  dr  KgV pg  pv  g 


W,  r  =  rx:  u  =  0,  v  =  0,  -Dh*  |ft-^-  =  0,  -Omr|,—  -0 


(26) 

(27) 

(28) 


(M- 


r  =  r2,  z 
^  i  071 

-^-+p  (c  +c  c)«r 

""  I*  Qr  'g  pg  pv  >  & 

r  ,  r5C 

-Z)  — +  mC 

“Ha  q 


I  Ot 


A  V 

^  ^  ,  3C 

-Z)J  — +  mC 

mr|o(  0r 


+  h, 


Jo, 


A'-AA-TA  t29) 


Jb  V 


'N  r)T 

+  Zm2(l-e6)(C*-C),  J  =  0,  -z>*  * 


dT 

-DtS-^  =  hs2(Tg-Ts) 
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z  =  0,  re(r2,r3):  u  =  0,  v  =  0,  -D,z8 


8T 

dC 

L t 

g 

=  o, 

-D\ 

—  =  0 

dz 

mz  \ 

ot  dz 

dT 

,  dC 

’ 

g 

=  0, 

—  Dm. 

\  —  =  0 

\Oi 

'  dz 

mi 

'  ot  dz 

(30) 

(31) 


Specific  boundary  conditions: 


For  the  1 -cartridge  adsorber,  the  basic  difference  between  the  three  models  is  in  the  boundary 
conditions  at  the  adsorber  wall  (for  r=r4): 

-  For  Model_A_l  and  Model_D_l  (adsorption  and  desorption  without  condensation)  they 
are  defined  as:  zero  velocity  at  the  wall,  zero  adsorbate  flux  through  the  wall  and  the  heat  flux  through 
the  wall  defined  by  the  heat  losses: 


dT 

—  =  0  - Dhg I  — -  =  h  (T  -T  ) 

-™\ot  ^  U’  Utr  \ot  ^  "wgVg  La) 


r  =  r4,  z  e.{h,h  +  H)\  u  =  0,  v  =  0,  -Dmr — =  0,  -D^\  ^  =  hwe(Tg-Ta)  (32a) 

-  For  Model  DC  l  (desorption  with  condensation)  the  boundary  conditions  at  the  wall 


become: 


r  =  r„  z  e(h,h  +  H):  u -0,  v  =  0,  C  =  C„(T„), -D* 


ot  dr 


,  dC 

-Dmr\ot  -^(.-^cond)  +  Kg(Tg -Ta){ 32b) 


meaning  that,  when  the  conditions  for  condensation  at  the  wall  are  fulfilled,  the  concentration  of  the 
gas  phase  at  the  wall  becomes  equal  to  the  saturation  concentration  of  the  vapor  corresponding  to  the 
wall  temperature  Tw,  Csat(Tw ).  At  the  same  time,  in  the  boundary  condition  for  the  heat  balance,  the 
heat  of  condensation  term  is  be  included  ((-A Hcond)  is  the  molar  heat  of  condensation).  In  our  model, 
the  heat  resistance  and  heat  capacity  of  the  wall  are  neglected  and  the  wall  temperature  is  assumed  to 
be  equal  to  the  ambient  temperature  Ta. 

B2)  Common  boundary  conditions  for  all  models  for  the  2-cartridges  adsorber 


z  =  0,  re(0,/i):  u  =  0,  v  =  — - ,  Tg=Tgin,  C  =  Cin 

'i  *P, 

z  =  hl+H,  re(0,r,):  u  =  0,  v  =  0,  -D*  1^  =  0,  -Dm  1^  =  0 


dz 


dz 


z  =  hl+  H  +  h2,  r  e  (0,r,)  :  u  —  0,  v  =  0,  -D?/\8f  =  0,  -Dmz\ff  =  0 

dz  u  dz 

z  g  (0,  ht ),  r  =  r,  :  u  =  0,  v  =  0,  -Dbg  |.(^-  =  0,  -Dmr \itf^  =  ° 


zt(hx,H\  r  =  rl:  p\it=p\b, 

(  ,  dT  ^  d 

■DhA.^y  +  Pg{cpg+cpvC)uTg 


u\u=u\b’  v\it=v\b’ 


\u  Qr 


dT 


g 


tr  \b  .  +P g(CPg+CpvC>Tg 
Jit  V  ur  Jb 


(33) 

(34) 

(35) 

(36) 


+  hsl(l-sb)(Ts-Tg) 


D  I  —  +  uC 

mrlit  dr 


Ju 


-Dmr\  —  +  uC 
mrlb  dr 


+  kml(\-&b)(C*  -C\ 


(37) 


Jb 


-Disfy=hsx(Tg-Ts\  j= o 

dT  dC  dT 

z  =  K,  r  e  (rl5r2)  :  u  =  0,  v  =  0, -D*\  =  0, -Dmz\  —  =  0,  U  =  U0, -Dbs  -*■  =  0  (38) 

16  dz  b  dz  dz 
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ST  SC  ST 

z  =  hl+H,  re(r„r2):  u  =  0,  v  =  0, -D^  -£■  =  0, -Dmz\b  —  =  0,  U  =  0, -D%  =  0  (39) 


ze(H  +  hl+h2,2H  +  hi+h2),  r  =  rl:  p\it=p\b,  «|tf=«|6,  v|,y=vL’ 

Z'  S77  ^  f  0J7  'N 

-^L-^  +  Pg(^+VC)M7;  =  -D?f\-^  +  Pg(cpg+cpvC)uTg  +hsl(l-sb)(Ts-Tg) 

V  1  ar  Jit  V  )b 

f  sc  \  f  sc  \ 

-Dmr\—  +  uC  =  -Dmr\—  +  uC  +  kmX(\-&b){C  -C), 

V  or  A-,  V  Or  Ju 


hJTg~Ts),  J  =  0 


ST  SC  ST 

z  =  hl+h2+H,  re(r„r2):  u  =  0,  v  =  0, -Z)*g|  -*■  =  0, —  =  0,  <7  =  <70, -D?  =  0(41) 

^  Sz  D  Sz  Sz 


ST  SC  ST 

z  -  2H  +  hx+h2,  re  (r„r2)  :  u  =  0,  v  =  0, -DhA  =  0,  -Z)mz|fi  —  =  0,  U  =  0,  -D%  =  0  (42) 

^  dz  0  Sz  Sz 

ze(OA),  r  =  r2:  u  =  0,  v  =  0,  -0*1  ^-  =  0,  -Dmr\—  =  0  (43) 


z  s(hl,hl+  H),  r  =  r2:  p\b=p\ot,  u\b=u\ot,  v|  =v|  , 


9T7  S  3T7  ^ 

~D‘r8l~dS' +P^Cps+CpT>TS  =\-D»l  -^r  +  PgOPg+cpT)uTg  +hs2(l-£b)(Ts  -Tg\ 

Jb  V  A;/ 


“  A»r  /,  “T“  +  “C  -  “A 
dr  A  V 


ST 


+  uC  +  km2(l  -  eb)(C*  -  C),  J  =  0, 


ze  (hx+H,hx+h2+H),  r  =  r2:  u=  0,  v  =  0,  -0^^  =  °’  -Dmr  |^  =  0 

zs(hl  +  h2+H,hl+h2  +  2H),  r  =  r2\  p\b=p\ot,  u\b=u\gt,  v|6=v|o(, 

,  dT  \  (  i  dT  \ 

~  D*l  if + p,(<W + ‘/H  =  h  °?L,  if + p,  <cK + +  Mi  -  s»xr,  -  rp, 

V  Jb  \  Jot 


mnb  dr 


~  Dmr\h^r  +  UC  \  ~  \~  Dmr  L  “Z-  +  +  “  8fe)(C  “  Q>  /  -  0, 


""•lot  Qr 


ST 

-D?-£r  =  h*V*-TJ 


z  =  2H  +  hx  +  h2,  r  e  (r2,r4) :  u  =  0,  v  =  0,  -0*g  — =  ° 

ot  dz  °‘  dz 


ST  ST  SC  SC 

z  e  (0,  hx),  re(r3,r4):  M  =  0,  v  =  0, -Z)*g|  — L-Dp\  — -  =  0,-Z>mf.|  — - Dmz\  — 

lo*  qz  i ot  sr  1 ot  Sr  ot  Sz 


■  0  (48) 


Specific  boundary  conditions: 

For  the  2-caartridges  adsorber,  the  differences  between  the  models  corresponding  to  adsorption  and 
desorption  are  a  result  of  different  gas  flow  directions  in  these  two  steps.  Accordingly,  the  boundary 
conditions  for  the  adsorber  inlets  and  outlets  are  different.  On  the  other  hand  the  boundary  conditions 
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at  the  adsorber  wall  differ  for  the  models  without  and  with  condensation  (in  an  analogous  way  as  for 
the  1 -cartridge  models).  The  corresponding  boundary  conditions  are: 

-For Model  A  2: 


'2>'3 ) 


~)hg 

Jtz 


I 

'ot  dz 


-  For  Model_D_2  and  Model  DC  2: 

.  .  dT 


I °t  dz 


K 

II 

© 

II 

o 

(49a) 

1 

^b 

II 

O 

1 

ib 

II 

o 

(50a) 

DC_2: 

a 

II 

(49b) 

z  =  hx+h2+2H,  u-  0,  v  =  - 


T  -T .  C=C 

g  gin  in 


(50b) 


-  For  Model_A_2  and  Model_D_2: 
ze(hi,2H  +  hl+h2),  r  =  r4:  u  =  0,  v  =  0,  ~Dmr\gt 
-For Model  DC  2: 


—  =  0  -Dhg I  -£-  =  h  (T-T) 

dz  ’  tr  \ot  dr  wgl  8  a) 


z  e  {hx,2H  +  hx  +  h2),  r  =  r4:  u  =  0,  v  =  0,  C  =  Csat(Tw), 

dC 


,  ,  dT„ 

-D*\-±  =  -Dm 
\Qt  or 


dr 


(-A Hcond)  +  h(T-Ta) 


(51a) 


(51b) 


In  equations  (1-51),  G  is  the  molar  flow-rate  of  the  gas  inert,  U0  -  the  supply  voltage  and  8b  -  bed 
porosity.  The  subscript  in  denotes  the  inlet,  o  the  outlet,  it  the  inner  (central)  tube,  ot  the  outer 
(annular)  tube  and  a  the  ambient.  The  meanings  of  the  geometric  values  ( rx  to  r4,  H,  h\  and  h2 )  are 
defined  in  Figure  5.  A  complete  list  of  the  notations,  including  the  physical  parameters  and 
coefficients  used  in  the  model  equations,  can  be  found  in  the  nomenclature,  at  the  end  of  this  report. 

Although  the  heat  and  mass  transfer  coefficients,  dispersion  coefficients,  etc.,  can  be  defined  as 
functions  of  the  dependent  variables  (velocities,  concentrations,  temperatures),  in  our  current  models 
constant  values  were  used  in  order  to  reduce  the  convergence  problems.  The  adsorption  isotherm 
relation  ( q=f(C ,  Ts))  and  the  temperature  dependence  of  the  electric  resistivity  were  included  in  the 
models  (equations  (22)  and  (16),  respectively). 
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3.  Models  for  complete  TSA  cycles 

In  order  to  model  a  complete  TSA  cycle,  it  is  necessary  to  switch  from  one  Femlab  model  to  another, 
when  certain  conditions  are  fulfilled. 

In  order  to  reduce  the  time  needed  to  achieve  the  periodic  quasi-steady  state,  the  TSA  cycle  was 
started  with  the  desorption  step.  The  initial  state  of  the  system  corresponds  to  equilibrium  at  ambient 
temperature.  The  initial  concentration  in  the  gas  phase  is  equal  to  the  feed  concentration  during 
adsorption: 

/  <  0,  Vr,  Vz :  Tg  =  Ts  =  Tgm  =  Tgo  =  Ta  =  Tp,  C  =  Cin  =  C*  =  Q>(q)  =  Cp,  u  =  v  =  0,  p  =  pa 


Time=21  [s]  Surface:  Concentration  Subdomain  marker:  concentration 


condensation 

started 


r 


r  [cm] 


Three  switching  moments  are  defined: 

1.  Switching  from  Model_D  to  ModelDC. 

This  switch  is  performed  at  the  moment  when 
the  gas  concentration  at  the  adsorber  wall  (for 
r=r4 )  becomes  equal  to  the  saturation 
concentration  corresponding  to  the  wall 
temperature  Csah  resulting  with  start  of 
condensation.  The  end  of  condensation  can  be 
detected  from  the  change  of  the  sign  of  the 
adsorbate  flux  at  the  adsorber  wall 
(- Dmr\ot(dC / dr\=r4  becomes  <0).  In  Figure 

11  we  show  the  distribution  of  the  gas 
concentrations  in  subdomain  COl  at  the  first 
switching  moment,  for  the  1-cartidge  adsorber. 

It  can  be  seen  that  the  concentration  is  nearly 
uniform  at  the  wall  surface,  meaning  that  the 
concentration  practically  starts  on  the  whole 

surface  simultaneously.  A  similar  result  is  obtained  for  the  2-cartridges  adsorber. 

2.  Switching  from  Model  DC  to  Model_A. 

This  switch  defines  the  end  of  the  desorption  and 
the  start  of  the  adsorption  part  of  the  cycle.  It  is 
performed  when  the  maximal  temperature  in  the 
adsorbent  bed  reaches  certain  predefined  value 
Tsw,  which  should  not  be  exceeded.  The 
temperature  distribution  in  the  solid  phase  at  the 
moment  when  the  maximal  temperature  reaches 
473.15  K  is  shown  in  Figure  12,  for  the  1-cartridge 
adsorber.  It  can  be  seen  that  significant 
temperature  gradients  are  present  both  in  the  axial 
and  in  the  radial  directions  in  the  adsorbent  bed. 


iiri  =  0.0104251 


■0.0104 


0.0104 


•  -U.U1U4 

U  0.01 04 
0.0104 
0.0104 


Figure  11.  Maximal  and  minimal  concentration  at  the 
adsorber  wall  in  the  1st  switching  moment  (1-catridge) 
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3.  Switching  from  Model  A  to  Model _D. 

This  switch  corresponds  to  the  end  of  adsorption 
and  start  of  the  new  desorption  step.  It  is 
performed  automatically  when  the  outlet 
concentration  reaches  a  certain  predefined  breakthrough  value 
concentration  is  obtained  as  the  average  value  across  the  outlet: 

1 


r  [cm] 

Figure  12.  Maximal  and  minimal  temperature  in  the 
adsorbent  bed  in  the  2nd  switching  moment  (1-catridge) 

Cbreak-  In  principle,  the  outlet 
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z=0 


dr . 
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Time=2884  [s]  Surface:  Concentration  Subdomain  marker:  concentration*  10'4 

From  Figure  13,  which  shows  the  concentrations 
in  the  annular  tube  of  the  1-catridge  adsorber,  at 
the  moment  of  breakthrough,  it  can  be  seen  that 
the  concentration  is  practically  uniform  across  the 
outlet.  The  same  result  is  obtained  for  the  2- 
cartridges  adsorber. 

Although  in  principle  different  Femlab  models  can 
be  run  consecutively,  automatic  checking  whether 
certain  conditions  are  met  and  switching  from  one 
model  to  another  based  on  it,  is  not  possible  in 
Femlab.  For  that  reason,  building  of  the  integral 
models  for  complete  TSA  cycles  was  done  by 
using  combination  of  Femlab  and  Matlab. 

The  procedure  for  building  these  models  can 
shortly  be  given  by  the  following  steps: 

1.  Building  and  running  the  Femlab  models  for  separate  stages  of  the  TSA  cycle  (Model_A_l, 
Model  D  l  and  Model  DC  l  for  1 -cartridge  and  Model_A_2,  Model_D_2  and  Model_DC_2 
for  2-cartridges  adsorber).  The  same  application  modes  have  to  be  used  in  all  models.  The 
same  sequence  of  the  application  modes  and  of  the  subdomains  has  to  be  used  in  all  models 
corresponding  to  the  same  type  of  adsorber. 

2.  Saving  the  Femlab  models  as  Matlab  programs  (m-files). 

3.  Integrating  all  three  m-files  into  one  Matlab  program  (some  program  segments,  like  geometry 
definition  and  mesh  initialization,  are  needed  only  once). 

4.  Adding  the  necessary  Matlab  commands: 

a.  For  checking  whether  certain  conditions  are  met; 

b.  For  defining  the  new  initial  state  as  the  solution  of  the  previous  step; 

c.  For  switching  from  one  to  another  model; 

d.  For  repeating  the  complete  cycle  the  needed  number  of  times. 

5.  Plotting  the  results  in  Matlab  or  exporting  them  as  a  fern  structure  that  can  be  imported  into 
Femlab  for  postprocessing. 

The  structure  of  the  integral  TSA  models  is  given  in  Figure  14.  This  structure  is  common  for  both 
configurations  of  the  adsorber  (1 -cartridge  and  2-cartridges),  so  common  names  Model  D,  Model  DC 
and  Model  A  are  used. 
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Figure  13.  Maximal  and  minimal  concentration  at  the  outlet 
in  the  3rd  switching  moment  (1-catridge) 
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Figure  14.  Structure  of  the  complete  TSA  model 
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4.  Simulation  results 


The  Femlab  and  Matlab/Femlab  models  were  used  for  simulation  of  adsorption,  electrothermal 
desorption  without  and  with  condensation  and  of  complete  TSA  cycles. 

Simulation  of  adsorption  was  performed  for  initially  clean  adsorbent  (qp= 0)  at  ambient  temperature, 
and  for  constant  feed  concentration  and  flow-rate. 

Simulation  of  desorption  was  performed  for  initially  uniformly  saturated  adsorbent  at  ambient 
temperature,  for  constant  voltage  supply  and  constant  flow-rate  inert  feed  stream  (Cin= 0).  The 
simulation  results  presented  in  this  report  correspond  to  a  complete  desorption  (regeneration)  step, 
including  desorption  before  the  start  of  condensation  and  desorption  accompanied  with  condensation. 

The  simulations  were  performed  for  a  laboratory  scale  adsorber  of  the  following  dimensions: 

-  ri=0.95  cm 

-  r2=  1.5  cm 

-  r3= 2  cm 

-  r4= 3.5  cm 

-  H=  30  cm 

-  h\=2  cm 

-  h2=2  cm 

and  for  the  following  system: 

-  Adsorbate:  methyl  ethyl  ketone  (MEK) 

-  Adsorbent:  American  Kynol  ACC-5092-20  (activated  carbon  fiber  cloth  material) 

-  Carrier  gas:  nitrogen. 

The  bed  permeability  was  calculated  from  the  experimentally  measured  pressure  drop  of  the  ACFC 
bed  given  in  Ref.  [5].  The  numerical  values  of  most  of  the  other  model  parameters  are  also  based  on 
Ref.  [5].  A  complete  list  of  the  model  parameters  used  in  these  simulations  is  given  in  Table  2. 

The  simulations  were  performed  with  the  initial  meshes  generated  in  Femlab,  and  shown  in  Figures 
10a  and  10b. 

Femlab  offers  a  large  number  of  options  for  graphical  representation  of  the  simulated  results.  The 
following  results  were  chosen  for  presentation  in  this  report: 

Velocity  field  in  the  adsorber; 

Gas  and  solid  concentration  fields  in  the  adsorber; 

Gas  and  solid  temperature  fields  in  the  adsorber. 

These  results  are  shown  in  different  forms:  as  2-D  surface  plots  and  as  time,  radial  and  axial  profiles, 
in  1-D  diagrams. 

For  desorption,  the  electric  resistivity  field  is  also  shown. 
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Table  2.  The  model  parameters  used  for  simulation 


Parameter 

Notation  and  value  | 

Adsorption 

Desorption  | 

Adsorber  dimensions 

^=0.95  cm  r2=1.5  cm 

r3= 2  cm  r4= 3.5  cm 

H=  30  cm  h\=  h2= 2  cm 

Bed  porosity 

86=0.72 

Specific  surface  area 

a=  13.65  cm2/cm3 

Ambient  pressure 

Patm=  101325  Pa 

Inlet  adsorbate  concentration  in  the  gas  phase 

C/n=0.001  mol/mol 

Cin= 0  mol/mol 

Initial  adsorbate  concentration  in  the  gas  phase 

Cp=0  mol/mol 

C^=0.001  mol/mol 

Initial  adsorbate  concentration  in  the  solid  phase 

qp=0  mol/g 

^=0.00465  mol/g 

Saturation  concentration 

Csat= 0.01  mol/mol  | 

Electric  voltage  supply  U0 

U0= 0  V 

U0= 11V 

Inert  gas  flow-rate 

G=0.1  mol/s 

G=0.01  mol/s 

Switching  temperature  from  desorption  to  adsorption 

7/^=463.15  K 

Breakthrough  concentration  for  switching  from  adsorption  to  desorption 

Cbreak=®- 0005  mol/mol 

Inlet  gas  temperature 

7^=293.15  K 

Ambient  temperature 

7fl=293.15  K 

Wall  temperature 

7^=293.15  K 

Mass  transfer  coefficient  in  the  bed 

£m=0.0003  mol/(cm2s) 

Bed  to  central  tube  mass  transfer  coefficient 

km 7=0.0003  mol/(cm2s) 

Bed  to  annular  tube  mass  transfer  coefficient 

km2=0.0 003  mol/(cm2s) 

Heat  transfer  coefficient  in  the  bed 

*,,=0.0055  W/(cm2s) 

Bed  to  central  tube  heat  transfer  coefficient 

*s7=0.001  W/(cm2s) 

Bed  to  central  tube  heat  transfer  coefficient 

*S2=0.001  W/(cm2s) 

Heat  transfer  coefficient  defining  heat  losses 

*„,e=5xl0'5  W/(cm2K) 

Radial  dispersion  coefficient  of  the  gas  phase  in  the  central  and  annular  tubes 

AJ  ,t=AJ  ot=l  mol/(cms) 

Axial  dispersion  coefficient  of  the  gas  phase  in  the  central  and  annular  tubes 

f^mzl  it- f^mzl  ot— 1  mol/(cms) 

Radial  heat  diffusivity  of  the  gas  phase  in  the  central  and  annular  tubes 

Arh8l  it=Arhgl  ot=l  W/(Kcm) 

Axial  heat  diffusivity  of  the  gas  phase  in  the  central  and  annular  tubes 

Azhgl  it=Azhgl  ot=l  W/(Kcm) 

Radial  and  axial  heat  diffusivity  of  the  solid  phase 

Afe=l-2xl0"3  W/(cmK) 

Molar  heat  of  adsorption 

(-A77atfe)=66200  J/mol 

Molar  heat  of  condensation 

(-A//cW)=3 1230  J/mol 

Adsorbent  bed  density 

p6=0.221  g/cm3 

Adsorbate  density 

p^=0.81  g/cm3 

Adsorbate  molar  mass 

Ma=12A01  g/mol 

Inert  gas  molar  mass  (nitrogen) 

Mb= 28.02  g/mol 

Specific  heat  capacity  of  liquid  adsorbate 

co/=157.9  J/(molK) 

Specific  heat  capacity  of  the  inert  gas 

<V=29.13  J/(molK) 

Heat  capacity  of  the  solid  phase 

Cm= 0.71  J/(gK) 

Specific  heat  capacity  of  gaseous  adsorbate 

CpV=  100.9  J/(molK) 

Dynamic  viscosity  of  the  inert  gas 

p=  1.78x1  O'4  g/cm/s 

Inert  gas  thermal  conductivity 

*=2.6xl0‘4W/cm/K 

Electric  resistivity  at  referent  temperature  TR 

p0=0.202  Hem 

Temperature  coefficient  of  the  bed  electrical  resistivity 

*=-2xl  0'3  1/K 

Critical  temperature  of  the  inert  gas  (nitrogen) 

TcB=  126.2  K 

Critical  pressure  of  the  adsorbate 

pcA= 4.26xlO+6Pa 

Critical  pressure  of  the  inert  gas  (nitrogen) 

/>cS=3.39xlO+6  Pa 

Wagner  constants  for  the  adsorbate 

VPA=- 7.71476,  TP5=1.71061 
EFc^-3.6877,  F7d=-0.75169 

Total  volume  of  micropores  (D-R  equation) 

1P0=0.748  cm3/g 

Boltzmann’s  constant 

*b=1.38048xl0'23  J/K 

Adsorption  energy  of  the  adsorbate  (D-R  eq.) 

£=14.43x1 0+3  J/mol 
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Two  additional  parameters,  very  important  for  the  electrothermal  desorption  process,  were  calculated: 
the  mass  of  the  condensed  liquid  and  the  used  electric  energy. 


Calculation  of  the  mass  of  the  condensed  liquid: 

The  flux  of  the  condensed  liquid  is  calculated  from  the  gas  concentration  gradient  at  the  adsorber  wall: 

(52) 

r=r4 

This  flux  is  a  function  of  the  axial  position  and  time.  The  condensation  rate  is  calculated  by  integration 
of  the  flux  over  the  whole  surface  of  the  adsorber  wall,  and  is  a  function  of  time  only: 

Htot 

Condi*)  =  2m4  \jcond(z,t)dz  (53) 

0 

The  total  amount  of  the  condensed  liquid  is  obtained  by  integrating  the  condensation  rate  over  time: 

^ cond 

LCond  ~  \LCond(f)dt  (54) 

0 

Htot  denotes  the  total  height  of  the  adsorber  and  Tcond  the  total  time  of  condensation. 


JCond(zd)  =  -Dn 


I 

lot  Qr 


Calculation  of  the  used  electric  energy: 

It  is  assumed  that  the  used  electric  energy  is  equal  to  the  Joule  heat  produced  in  the  adsorbent  material. 
This  energy  per  unit  volume  and  unit  time  is  defined  by  equation  (18). 

The  used  electric  power  is  obtained  by  integration  over  the  adsorbent  bed  volume: 


Htot  r=r3  1  ( 

fr  1 

(  dU ) 

2 

+ 

f 

2  A 

Lf  P0(l  +  b(Ts-TR)){ 

l  dr  J 

[dz  j 

J 

2ndrdz 


(55) 


and  the  total  used  energy  during  desorption,  by  integration  over  time: 


T des 

Qel  =  \Qeli*)dt 
0 

x des  is  the  total  desorption  time. 


(56) 


4.1.  Simulation  results  for  the  1-cartridge  adsorber 

4.1.1 .  Adsorption  in  the  1-cartridge  adsorber 

The  simulation  of  adsorption  is  performed  for  7000  s.  After  that  time  the  column  is  practically 
saturated  and  the  adsorption  is  finished.  The  results  presented  here  were  obtained  for  inlet 
concentration  0.001  mol/mol  and  inert  gas  flow  rate  0.1  mol/s. 

Velocity  distribution 

The  velocity  field  in  the  1 -cartridge  adsorber  is  presented  first.  In  Figure  15,  the  velocity  distribution 
at  the  end  of  adsorption  (at  t= 7000  s)  is  shown  in  the  form  of  a  2-D  surface  plot,  together  with  a  vector 
representation  of  the  velocity.  Constant  velocity  field  was  obtained  nearly  instantaneously  (during  the 
first  3  seconds).  The  velocity  distribution  is  also  presented  in  1-D  diagrams:  the  time  profiles  at  the 
middle  cross-section  of  the  adsorber  (z=16  cm)  and  4  radial  positions  in  Figure  16,  the  radial  profiles 
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for  *=7000  s  and  4  axial  positions  in  Figure  17  and  the  axial  profiles  for  *=7000  s  and  4  radial  positions 
in  Figure  18.  Some  numerical  noise  is  observed  in  the  simulated  gas  velocities  in  the  central  tube. 


Time=7000  s  Surface:  wg 


r  [rml 


Figure  15.  Velocity  distribution  of  the  gas 
phase  during  adsorption  in  the  1  -cartridge 
adsorber,  for  £=7000s 


Figure  17.  Radial  profiles  of  the  gas  velocity  during 
adsorption  in  the  1 -cartridge  adsorber,  for  £=7000  s 


Figure  16.  Time  profiles  of  the  gas  velocity  during 
adsorption  in  the  1 -cartridge  adsorber,  at  z=  16  cm 


Figure  18.  Axial  profiles  of  the  gas  velocity  during 
adsorption  in  the  1 -cartridge  adsorber,  for  £=7000  s 


Concentration  distributions 

The  gas  and  solid  concentration  distributions  during  adsorption  in  the  1 -cartridge  adsorber  are  first 
presented  in  the  form  of  2-D  surface  plots,  in  Figures  19  and  20  respectively.  The  concentrations  at 
three  distinct  moments  (*=50,  200  and  7000  s)  are  shown,  in  order  to  illustrate  their  evolution  in  time. 
For  *=7000  s  the  gas  concentration  in  the  adsorber  is  nearly  uniform  and  equal  to  the  feed 
concentration  (0.001  mol/mol),  which  means  that  the  adsorption  process  is  finished.  At  early  stages  of 
adsorption,  the  concentrations  in  the  adsorber  are  nonuniform,  both  in  the  axial  and  in  the  radial 
direction. 

The  concentration  profiles  are  also  given  in  1-D  diagrams:  the  time  profiles  at  z=16  cm  in  Figure  21, 
the  radial  profiles  at  z=16  cm  in  Figure  22  and  the  axial  profiles  for  *=200s  in  Figure  23.  Figures  21a, 
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22a  and  23a  correspond  to  the  gas  phase  and  Figures  21b,  22b  and  23b  to  the  solid  phase.  For  the 
model  parameters  used  in  our  simulations,  the  gas  concentration  radial  profiles  in  the  central  and  in  the 
annular  tube  are  practically  flat. 


(a)  (b)  (c) 
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Figure  19.  Gas  concentration  distribution  during  adsorption  in  the  1 -cartridge  adsorber:  (a)  £=50s,  (b)  /=200s  and 
(c)  t=7000s 
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Figure  20.  Solid  concentration  distribution  during  adsorption  in  the  1 -cartridge  adsorber:  (a)  £=50s,  (b)  /=200s  and 
(c)  *=7000s 
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Figure  21.  Time  profiles  of  the  concentrations  during  adsorption  in  the  1 -cartridge  adsorber,  at  z=16cm 
(a)  gas  phase;  (b)  solid  phase 


hgure  22.  Radial  profiles  of  the  concentrations  during  adsorption  in  the  1 -cartridge  adsorber,  at  z=16cm: 
(a)  gas  phase;  (b)  solid  phase 


(a)  (b) 


Figure  23.  Axial  profiles  of  the  concentrations  during  adsorption  in  the  1 -cartridge  adsorber,  for  £=200  s: 
(a)  gas  phase;  (b)  solid  phase 
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Temperature  distributions 

The  gas  and  solid  temperature  distributions  in  the  1 -cartridge  adsorber  are  first  given  in  Figures  24  and 
25,  respectively,  in  the  form  of  2-D  surface  plots.  The  simulated  gas  and  solid  temperatures  are  also 
shown  in  1-D  diagrams:  the  time  profiles  in  Figures  26a  and  26b,  the  radial  profiles  in  Figures  27a  and 
27b  and  the  axial  profiles  in  Figures  28a  and  28b.  All  these  figures  show  temperature 
nonhomogeneity  in  both  phases.  As  the  temperature  changes  in  the  adsorber  are  due  only  to  the  heat 
generated  by  adsorption,  they  are  not  significant  (the  maximal  temperature  change  is  less  than  5  °C. 
The  biggest  temperature  change  is  observed  during  the  first  200  s  of  the  adsorption  process. 


(a) 

Time=50  s  Surface:  Gas  temperature 


(b) 

Time=200  s  Surface:  Gas  temperature 
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Figure  24.  Gas  temperature  distribution  during  adsorption  in  the  1 -cartridge  adsorber:  (a)  £=50s,  (b)  £=200s  and  (c)  £=7000s 
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fime=50  s  Surface:  Solid  temperature 
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Time=200  s  Surface:  Solid  temperature 
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Figure  25.  Solid  temperature  distribution  during  adsorption  in  the  1 -cartridge  adsorber:  (a)  t=50s,  (b)  t=200s  and  (c)  t=7000s 
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Figure  26.  Time  profiles  of  the  temperatures  during  adsorption  in  the  1 -cartridge  adsorber,  atz=16cm: 
(a)  gas  phase;  (b)  solid  phase 
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Figure  26.  Radial  profiles  of  the  temperatures  during  adsorption  in  the  1 -cartridge  adsorber,  atz=16cm: 
(a)  gas  phase;  (b)  solid  phase 
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Figure  28.  Axial  profiles  of  the  temperatures  during  adsorption  in  the  1 -cartridge  adsorber,  for  t= 200  s 
(a)  gas  phase;  (b)  solid  phase 
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4.1.2.  Desorption  in  the  1-cartridge  adsorber 

The  complete  desorption  step,  before  and  after  the  start  of  condensation,  was  simulated.  This  was  done 
by  using  a  Matlab  model  obtained  by  combining  Model_D_l  and  Model_DC_l.  The  results  presented 
here  were  obtained  for  initial  concentration  in  the  solid  phase  0.00465  mol/g  (equilibrium 
concentration  for  C=0.001  mol/mol),  electric  voltage  11  V  and  inert  gas  flow-rate  0.01  mol/s.  The 
time  scale  used  for  this  simulation  was  0  to  200  s.  After  200  s  the  desorption  process  is  nearly 
finished. 


Velocity  distribution 

The  velocity  distribution  at  the  end  of  desorption  (at  £=200  s)  is  shown  in  Figure  29,  in  the  form  of  a  2- 
D  surface  plot,  together  with  a  vector  representation  of  the  velocity.  The  velocity  distribution  is  also 
presented  in  1-D  diagrams:  the  time  profiles  in  Figure  30,  the  radial  profiles  in  Figure  31  and  the  axial 
profiles  in  Figure  32.  Some  (not  significant)  change  of  the  velocity  profile  with  time  is  observed,  as  a 
result  of  the  change  of  the  gas  temperature. 


Time=200s  Surface:  Gas  velocity 
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Figure  29.  Velocity  distribution 
during  desorption  in  the  1 -cartridge 
adsorber,  for  t=200  s 


Figure  30.  Time  profiles  of  the  gas  velocity  during 
desorption  in  the  1 -cartridge  adsorber,  at  z=  16  cm 


Figure  31.  Radial  profiles  of  the  gas  velocity  during 
desorption  in  the  1 -cartridge  adsorber,  for  t=200  s 


Figure  32.  Axial  profiles  of  the  gas  velocity  during 
desorption  in  the  1 -cartridge  adsorber,  for  t=  200  s 
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Concentration  distributions 

The  gas  and  solid  concentration  distributions  during  desorption  in  the  1 -cartridge  adsorber,  presented 
in  the  form  of  2-D  surface  plots,  are  given  in  Figures  33  and  34  respectively.  The  concentrations  at 
three  distinct  moments  (£=24,  80  and  200  s)  are  shown.  £=24  s  is  the  moment  when  condensation 
starts,  i.e.  when  the  program  switches  from  Model  D  l  to  Model  DC  l.  After  that  time  the 
concentration  at  the  adsorber  wall  (and  practically  in  the  whole  annular  tube)  becomes  constant  and 
equal  to  the  saturation  concentration  at  293.15  K  (0.01  mol/mol  for  MEK). 

The  concentration  profiles  are  also  given  in  1-D  diagrams:  the  time  profiles  at  z=16  cm  in  Figure  35, 
the  radial  profiles  at  z=16  cm  in  Figure  36  and  the  axial  profiles  for  t  =  80  s  in  Figure  37.  Figures  35a, 
36a  and  37a  correspond  to  the  gas  phase  and  Figures  35b,  36b  and  37b  to  the  solid  phase. 
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Figure  33.  Gas  concentration  distribution  during  desorption  in  the  1 -cartridge  adsorber:  (a)  t=  24s,  (b)  t=  80s  and 
(c)  £=200s 
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Figure  34.  Solid  concentration  distribution  during  desorption  in  the  1 -cartridge  adsorber:  (a)  t= 24s,  (b)  £=80s  and 
(c)  /=200s 
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Figure  35.  Time  profiles  of  the  concentrations  during  desorption  in  the  1-cartridge  adsorber,  atz=16cm: 
(a)  gas  phase;  (b)  solid  phase 
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Figure  36.  Radial  profiles  of  the  concentration  during  desorption  in  the  1 -cartridge  adsorber,  atz=16cm: 
(a)  gas  phase;  (b)  solid  phase 


(a)  (b) 

x  10'3 


Figure  37.  Axial  profiles  of  the  concentrations  during  desorption  in  the  1 -cartridge  adsorber,  for  /=80  s: 
(a)  gas  phase;  (b)  solid  phase 
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Condensation  flux  and  condensation  rate 

The  condensation  flux  is  calculated  using  equation  (52).  It  is  a  function  of  both  time  and  axial  position 
in  the  adsorber.  The  time  profiles  of  the  condensation  flux  are  shown  in  Figure  38,  and  the  axial 
profiles  in  Figure  39. 


Figure  38.  Time  profiles  of  the  condensation  flux  in 
the  1  -cartridge  adsorber 


These  figures  show  that  condensation  starts  at 
t= 24  s  and  the  condensation  flux  is  maximal  at 
approximately  £=80s.  At  £=200s,  when  the  power 
is  switched  of,  the  condensation  is  still  not 
completely  finished.  At  the  beginning  of  the 
process  the  condensation  flux  increases  from  the 
bottom  to  the  top  of  the  adsorber,  but  with  time 
the  axial  profile  changes  and  the  condensation  is 
most  intensive  near  the  bottom  of  the  adsorber. 

Figure  40  shows  the  time  profile  of  the 
condensation  rate,  calculated  by  integration  of 
the  condensation  flux  along  the  adsorber  wall 
surface  (equation  (55)). 


Figure  39.  Axial  profiles  of  the  condensation  flux  in 
the  1  -cartridge  adsorber 


TIME  (s) 


Figure  40.  Condensation  rate  in  the  1 -cartridge  adsorber 


Temperature  distributions 

The  2_D  surface  plots  representing  the  gas  and  solid  temperature  distributions  in  the  1 -cartridge 
adsorber  are  given  in  Figures  41  and  42,  respectively.  The  simulated  gas  and  solid  temperatures  are 
also  shown  in  1-D  diagrams:  the  time  profiles  in  Figures  43a  and  43b,  the  radial  profiles  in  Figures 
44a  and  44b  and  the  axial  profiles  in  Figures  45a  and  45b.  All  these  figures  show  temperature 
nonhomogeneity  in  both  phases  and  both  directions.  During  desorption,  the  temperatures  in  the 
adsorber  change  significantly,  owing  to  the  production  of  Joule  heat. 
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Figure  41.  Gas  temperature  distribution  during  desorption  in  the  1 -cartridge  adsorber:  (a)  J=24s,  (b)  J=80s  and 
(c)  /=200s 
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Figure  42.  Solid  temperature  distribution  during  desorption  in  the  1 -cartridge  adsorber:  (a)  £=24s,  (b)  J=80s  and  (c)  £=200s 
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Figure  43.  Time  profiles  of  the  temperatures  during  desorption  in  the  1-cartridge  adsorber,  atz=16cm: 
(a)  gas  phase;  (b)  solid  phase 
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Figure  44.  Radial  profiles  of  the  temperatures  during  desorption  in  the  1 -cartridge  adsorber,  at  z=16cm: 
(a)  gas  phase;  (b)  solid  phase 


(a)  (b) 


Figure  45.  Axial  profiles  of  the  temperatures  during  desorption  in  the  1 -cartridge  adsorber,  at  r=1.25  cm: 
(a)  gas  phase;  (b)  solid  phase 
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Distribution  of  the  electric  resistivity  of  the  ACFC  adsorbent  bed 

The  electric  resistivity  of  the  ACFC  is  temperature  dependent  (decreases  with  the  increase  of 
temperature).  As  the  temperature  in  the  adsorbent  bed  varies  considerably  both  in  time  and  space,  so 
does  the  electric  resistivity.  This  dependence  is  shown  in  Figure  46  (in  the  form  of  a  2-D  surface  plot) 
and  Figures  47,  48  and  49  (1-D  diagrams  showing  time,  radial  and  axial  distributions,  respectively). 

The  electric  power  transformed  into  Joule  heat  also  changes,  as  a  consequence  of  the  change  of  the 
electric  resistance.  For  constant  voltage  supply,  the  electric  power  increases  with  time,  as  shown  in 
Figure  50.  The  electric  power  is  calculated  using  equation  (55). 
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Figure  46.  The  electric  resistivity  of  the  ACFC  adsorbent  distribution  in  the  1 -cartridge  adsorber:  (a)  £=24s,  (b)  *=80s  and 
(c)  /=200s 


Figure  48.  Radial  profiles  of  the  electric  resistivity 

Figure  47.  Time  profiles  of  the  electric  resistivity  during  desorption  in  the  1 -cartridge  adsorber,  at 

during  desorption  in  the  1 -cartridge  adsorber,  at  z=16cm 

z=  1 6cm 
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Figure  49.  Axial  profiles  of  the  electric  resistivity  in 
the  1 -cartridge  adsorber,  for  £=80  s 


Figure  50.  Electric  power  supply  during  desorption 
in  the  1  -cartridge  adsorber 


4.1.3.  TSA  process  in  the  1-cartridge  adsorber 

The  complete  TSA  process  was  simulated  by 
using  the  integral  Matlab/Femlab  model  presented 
in  Figure  14.  These  simulations  take  long  times,  so 
they  were  limited  to  the  first  three  cycles.  The 
simulation  results  presented  here  were  obtained  for 
r5W=463.15  K  and  Cbreak=  0.5  Cin.  All  other 
simulation  parameters  are  the  same  as  in  the 
simulations  of  adsorption  and  desorption.  The  time 
and  space  distributions  of  the  velocity, 
concentrations  and  temperatures  in  the  system 
during  the  adsorption  and  desorption  half-cycles 
are  very  similar  to  the  results  presented  in  Sections 
4.1.1  and  4.1.2,  so  they  will  not  be  given  here. 
Instead,  we  are  presenting  the  change  of  the  outlet 
concentration  from  the  adsorber  during  the  first  3 
cycles  in  Figure  51,  the  maximal  temperature  in  the 
rate  in  Figure  53.  These  figures  show  that  the  qu 
second  cycle. 


Figure  51.  The  outlet  concentration  during  TSA  in 
the  1 -cartridge  adsorber  (3  cycles) 


adsorbent  bed  in  Figure  52  and  the  condensation 
asi-steady  state  is  practically  achieved  after  the 


Figure  52.  The  maximal  adsorbent  temperature 
during  TSA  in  the  1 -cartridge  adsorber  (3  cycles) 


Figure  53.  The  condensation  rate  during  TSA  in  the  1- 
cartridge  adsorber  (3  cycles) 
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In  Figure  53  an  enlarged  segment,  showing  the  desorption  step  during  the  second  cycle,  is  also  given. 
From  this  segment,  it  can  be  seen  that  at  the  end  of  the  desorption  step,  the  condensation  hasn’t  been 
finished. 


4.2.  Simulation  results  for  the  2-cartridges  adsorber 
4.2.1.  Adsorption  in  the  2-cartridges  adsorber 

Owing  to  larger  amount  of  adsorbent,  the  time  needed  to  complete  the  adsorption  process  in  the  2- 
cartridges  adsorber  is  much  longer  than  for  the  1 -cartridge  adsorber.  For  that  reason,  these  simulations 
were  run  for  11000  s  or  ~  3  hours.  The  results  presented  here  were  obtained  for  the  same  inlet 
concentration  and  flow-rate  as  for  the  1 -cartridge  adsorber  (0.001  mol/mol  and  0.1  mol/s). 

Velocity  distribution 

The  velocity  field  in  the  2-cartridges  adsorber  is  presented  in  Figures  54-57.  In  Figure  54,  the  velocity 
distribution  at  the  end  of  adsorption  (at  t=l  10000  s)  is  shown  in  the  form  of  a  2-D  surface  plot, 
together  with  a  vector  representation  of  the  velocity.  The  velocity  distribution  is  also  presented  in  1-D 
diagrams:  the  time  profiles  at  the  middle  cross-sections  of  both  cartridges  (z=16  cm  andz=42  cm)  and 
4  radial  positions  in  Figure  55,  the  radial  profiles  at  z=16  cm  and  z=42  cm  in  Figure  56  and  the  axial 
profiles  in  several  radial  positions  in  Figure  57. 


Time=1 1000  s  Surface:  wg 
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Figure  54.  Velocity  distribution  of  the 
gas  phase  during  adsorption  in  the  2- 
cartridge  adsorber  for  t= 1 1000s 


Figure  55.  Time  profiles  of  the  gas  velocity  during 
adsorption  in  the  2-cartridge  adsorber  at  z=  16  cm 
(solid  lines)  and  z=49  cm  (dashed  lines) 
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Figure  56.  Radial  profiles  of  the  gas  velocity 
during  adsorption  in  the  2-cartridge  adsorber  for 
t=  11000  s 


Figure  57.  Axial  profiles  of  the  gas  velocity  during 
adsorption  in  the  2-cartridge  adsorber  at  r=1.25  cm 


Concentration  distributions 

The  2-D  surface  plots  of  the  gas  and  solid  concentration  distributions  during  adsorption  in  the  2- 
cartridges  adsorber  are  given  in  Figures  58  and  59,  respectively.  The  concentrations  at  three  distinct 
moments  (£=100,  400  and  11000  s)  are  shown.  For  £=11000  s  the  gas  concentration  in  the  adsorber  is 
uniform  and  equal  to  the  feed  concentration  (0.001  mol/mol),  i.e.,  the  adsorption  process  is  finished. 
At  early  stages  of  adsorption,  the  concentrations  in  the  adsorber  are  nonuniform. 

The  concentration  profiles  are  also  given  in  1-D  diagrams:  the  time  profiles  at  z=16  cm  in  Figure  60, 
the  radial  profiles  at  z=16  cm  in  Figure  61  and  the  axial  profiles  for  £=200  s  in  Figure  62.  Figures  60a, 
61a  and  62a  correspond  to  the  gas  phase  and  Figures  60b,  61b  and  62b  to  the  solid  phase. 
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Time:100  s  Surface:  Gas  concentration  Time=400  s  Surface:  Gas  concentration  Time-11000  s  Surface:  Gas  concent  y 


Figure  60.  Gas  concentration  distribution  during  adsorption  in  the  2-cartridges  adsorber:  (a)  £=100s,  (b)  /=400s  and 
(c)  t=  \  1000s 
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Figure  61.  Solid  concentration  distribution  during  adsorption  in  the  2-cartridges  adsorber:  (a)  £=100s,  (b)  /=400s  and 
(c)*=l  1000s 


Figure  62.  Time  profiles  of  the  concentrations  in  the  2-cartridges  adsorber,  at  z=16cm  (solid  lines)  and  at  z=49cm 
(dashed  lines):  (a)  gas  phase;  (b)  solid  phase 
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Figure  63.  Radial  profiles  of  the  concentrations  during  adsorption  in  the  2-cartridges  adsorber,  atz=16cm  (solid 
lines)  and  z=49cm  (dashed  lines):  (a)  gas  phase;  (b)  solid  phase 


38 


Figure  64.  Axial  profiles  of  the  concentrations  during  adsorption  in  the  2-cartridges  adsorber,  for  t= 200  s: 

(a)  gas  phase;  (b)  solid  phase 

Considerable  differences  between  the  concentrations  in  the  two  cartridges  can  be  observed.  These 
differences  diminish  in  time,  as  both  cartridges  become  nearly  saturated. 


Temperature  distributions 

The  gas  and  solid  temperature  distributions  during  adsorption  in  the  2-cartridges  adsorber  are  given  in 
Figures  65  and  66,  respectively,  in  the  form  of  2-D  surface  plots.  The  simulated  temperatures  are  also 
shown  in  1-D  diagrams:  the  time  profiles  in  Figures  67a  and  67b,  the  radial  profiles  in  Figures  68a  and 
68b  and  the  axial  profiles  in  Figures  69a  and  69b.  All  these  figures  show  temperature 
nonhomogeneity  in  both  phases.  The  maximal  temperature  rise  during  adsorption  is  less  than  4  °C,  for 
the  2-cartridges  adsorber.  The  highest  temperature  is  observed  after  approximately  400  s  from  the  start 
of  adsorption. 
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Figure  65.  Gas  temperature  distribution  during  adsorption  in  the  2-cartridges  adsorber:  (a)  £=100s,  (b)  /=400s  and 
(c)/=l  1000s 
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Figure  66.  Solid  temperature  distribution  during  adsorption  in  the  2-cartridges  adsorber:  (a)  t=l00s,  (b)  /=400s  and 
(c)*=l  1000s 
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Figure  67.  Time  profiles  of  the  temperatures  during  adsorption  in  the  2-cartridges  adsorber,  at  z=16cm  (solid  lines) 
and  z=49cm  (dashed  lines):  (a)  gas  phase;  (b)  solid  phase 


(a)  (b) 


Figure  68.  Radial  profiles  of  the  temperatures  during  adsorption  in  the  2-cartridges  adsorber,  at  z=16cm  (solid 
lines)  and  z=49cm  (dashed  lines):  (a)  gas  phase;  (b)  solid  phase 
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Figure  69.  Axial  profiles  of  the  temperatures  during  adsorption  in  the  2-cartridges  adsorber,  for  t= 200  s: 
(a)  gas  phase;  (b)  solid  phase 

Large  differences  between  the  temperature  profiles  in  the  two  cartridges  can  also  be  observed. 


4.2.2.  Desorption  in  the  2-cartridges  adsorber 

Again,  the  complete  desorption  step  is  simulated,  by  using  a  combination  of  Model_D_2  and 
Model_DC_2.  The  results  presented  here  were  obtained  for  initial  concentration  in  the  solid  phase 
0.00465  mol/g,  electric  voltage  11  V  and  inert  gas  flow-rate  0.01  mol/s  per  cartridge  (0.002  mol/s  in 
total).  The  time  scale  used  for  this  simulation  is  the  same  as  for  the  1 -cartridge  adsorber:  0  to  200s. 


Velocity  distribution 

The  velocity  distribution  at  the  end  of  desorption  (at  £=200  s)  is  shown  in  Figure  70,  in  the  form  of  a  2- 
D  surface  plot,  together  with  vector  representation  of  the  velocity.  The  velocity  distribution  is  also 
presented  in  1-D  diagrams:  the  time  profiles  in  Figure  71,  the  radial  profiles  in  Figure  72  and  the  axial 
profiles  in  Figure  73.  Some  (not  significant)  change  of  the  velocity  profile  in  time  is  present,  as  a 
result  of  the  change  of  the  gas  temperature. 
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Figure  70.  Velocity  distribution  of  the 
gas  phase  during  desorption  in  the  2- 
cartridge  adsorber  at  £=200  s 


Figure  71.  Time  profiles  of  the  gas  velocity  during 
desorption  in  the  2-cartridge  adsorber,  at  z=  16cm 
(solid  lines)  and  at  z=49cm  (dashed  lines) 
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Figure  72.  Radial  profiles  of  the  gas  velocity  Figure  73.  Axial  profiles  of  the  gas  velocity  during 

during  desorption  in  the  2-cartridge  adsorber,  for  desorption  in  the  2-cartridge  adsorber,  for  t=  200  s 

t=200  s 


Concentration  distributions 

The  gas  and  solid  concentration  distributions  during  desorption  in  the  2-cartridges  adsorber,  presented 
in  the  form  of  2-D  surface  plots,  are  given  in  Figures  74  and  75,  respectively.  The  concentrations  at 
three  distinct  moments  (7=30,  90  and  200  s)  are  shown  (7=30  s  is  the  moment  when  condensation 
starts,  i.e.  when  the  program  switches  from  Model_D_2  to  Model_DC_2).  The  concentration  profiles 
are  also  shown  in  1-D  diagrams:  the  time  profiles  in  Figure  76,  the  radial  profiles  in  Figure  77  and  the 
axial  profiles  in  Figure  78.  Figures  76a  to  78a  correspond  to  the  gas  phase  and  Figures  76b  to  77b  to 
the  solid  phase. 
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(c) 


Time=30s  Surface:  Gas  concentration 


Time=90s  Surface:Gas  concentration 
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Figure  74.  Gas  concentration  distribution  during  desorption  in  the  2-cartridges  adsorber:  (a)  /=30s,  (b)  /=90s  and  (c) 
/=200s 


Contrary  to  the  case  of  adsorption,  for  desorption  the  concentration  distributions  in  the  two  cartridges 
are  similar.  This  is  a  result  of  the  symmetrical  inlet  streams  in  the  two  cartridges. 
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(a)  (b)  (c) 

Time=30s  Surface:  Solid  concentration  Time=90s  Surface:  Solid  concentration  Time=200s  Surface:  Solid  concentration  10'3 


15 

Figure  75.  Solid  concentration  distribution  during  desorption  in  the  2-cartridges  adsorber:  (a)  £=30s,  (b)  *=90s  and  (c)  /=200s 
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Figure  76.  Time  profiles  of  the  concentrations  during  desorption  in  the  2-cartridge  adsorber,  atz=16cm  (solid 
lines)  and  at  z=49cm  (dashed  lines):  (a)  gas  phase;  (b)  solid  phase 


(a)  (b) 


Figure  77.  Radial  profiles  of  the  concentrations  during  desorption  in  the  2-cartridge  adsorber,  atz=16cm  (solid 
lines)  and  at  z=49cm  (dashed  lines):  (a)  gas  phase;  (b)  solid  phase 
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Figure  78.  Axial  profiles  of  the  concentrations  during  desorption  in  the  2-cartridge  adsorber,  for  J=80  s: 
(a)  gas  phase;  (b)  solid  phase 


Condensation  flux  and  condensation  rate 

The  condensation  flux,  calculated  by  using  equation  (52)  is  shown  in  Figure  79  (the  time  profiles)  and 
Figure  80  (the  axial  profiles). 


x10‘7 


Figure  79.  Time  profiles  of  the  condensation  flux  in 
the  2-cartridge  adsorber 


x  io'3 


Figure  81.  Condensation  rate  in  the  2-cartridges  adsorber 


Figure  80.  Axial  profiles  of  the  condensation  flux  in 
the  2-cartridges  adsorber 


These  figures  show  that  condensation  starts  at 
*=30  s  and  the  condensation  flux  is  maximal  at 
approximately  *=60s.  At  *=200s,  when  the 
power  is  switched  of,  the  condensation  is  still 
not  completely  finished.  At  the  beginning  of  the 
process  the  condensation  flux  increases  from 
the  adsorber  ends  towards  its  middle,  with  a 
sharp  minimum  at  the  middle  cross-section, 
between  the  two-cartridges. 

Figure  81  shows  the  time  profile  of  the 
condensation  rate,  calculated  by  integration  of 
the  condensation  flux  along  the  adsorber  wall 
surface. 
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Temperature  distributions 

The  gas  and  solid  temperature  distributions  in  the  2-cartridges  adsorber  in  which  desorption  is  taking 
place  are  again  presented  in  the  form  of  2-D  surface  plots  (Figures  82  and  83)  and  in  1-D  plots:  the 
time  profiles  in  Figures  84a  and  84b,  the  radial  profiles  in  Figures  85a  and  85b  and  the  axial  profiles  in 
Figures  86a  and  86b.  All  these  figures  show  temperature  nonhomogeneity  in  both  phases  and  large 
temperature  changes  during  desorption.  Similarity  between  the  temperature  distributions  in  the  two 
cartridges  can  easily  be  observed. 


(a) 


(b) 


(c) 


Time=30s  Surface:  Gas  temperature  Time=90s  Surface:  Gas  temperature  Time=200s  Surface:  Gas  temperature 
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Figure  82.  Gas  temperature  distribution  during  desorption  in  the  2-cartridges  adsorber:  (a)  *=30s,  (b)  *=90s  and  (c)  t=200s 


w  (b)  (c) 

Time=30s  Surface:  Solid  temperature  Time=90s  Surface:  Solid  temperature  Time=200s  Surface:  Solid  temperature 
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Figure  83.  Solid  temperature  distribution  during  desorption  in  the  2-cartridges  adsorber:  (a)  /=30s,  (b)  t= 90s  and  (c)  t= 200s 
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Gas  temperature  [K]  Gas  temperature  [K]  Gas  temperature  [K] 
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(b) 


Figure  84.  Time  profiles  of  the  temperatures  during  desorption  in  the  2-cartridges  adsorber,  atz=16cm  (solid  lines) 
and  z=49cm  (dashed  lines):  (a)  gas  phase;  (b)  solid  phase 


(a)  (b) 


Figure  85.  Radial  profiles  of  the  temperatures  during  desorption  in  the  2-cartridges  adsorber,  at  z=16cm  (solid 
lines)  and  z=49cm  (dashed  lines):  (a)  gas  phase;  (b)  solid  phase 


(a)  (b) 


Figure  86.  Axial  profiles  of  the  temperatures  during  desorption  in  the  2-cartridges  adsorber,  for  J=80  s 
(a)  gas  phase;  (b)  solid  phase 
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Electric  resistivity  [Qcm] 


Distribution  of  the  electric  resistivity  of  the  ACFC  adsorbent  bed 

The  distribution  of  the  electric  resistivity  in  the  2-cartridges  adsorber  is  shown  in  Figure  87  (in  the 
form  of  2-D  surface  plot)  and  in  Figures  88,  89  and  90  (the  1-D  diagrams  showing  the  time,  radial  and 
axial  distributions,  respectively). 
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Figure  87.  Distribution  of  the  electric  resistivity  of  the  ACFC  adsorbent  in  the  2-cartridges  adsorber:  (a)  /=30s,  (b)  /=90s  and 
(c)  /=200s 


Figure  88.  Time  profiles  of  the  electric  resistivity  in 
the  2-cartridges  adsorber,  at  z=  16cm  (solid  lines)  and 
at  z=49cm  (dashed  lines') 


Figure  89.  Radial  profiles  of  the  electric  resistivity  in 
the  2-cartridges  adsorber,  at  z=  16cm  (solid  lines)  and 
z=49cm  (dashed  lines') 


The  time  dependence  of  the  electric  power  is  shown  in  Figure  91. 
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Maxi  mal  adsorbent  temperatu  re  [K]  Outlet  concentrtion  [mol/mol] 


Figure  90.  Axial  profiles  of  the  electric  resistivity  in 
the  2-cartridges  adsorber,  for  J=80  s 


4.2.3.  TSA  process  in  the  2-cartridges  adsorber 


Figure  92.  Time  profiles  of  the  outlet  concentration 
in  the  2-cartridges  adsorber,  for  3  cycles 


Figure  93.  Maximal  adsorbent  temperature  in  the 
2-cartridges  adsorber  for  3  cycles 


Figure  91.  Used  electric  energy  during  desorption  in  the  2- 
cartridges  adsorber 


The  simulation  results  presented  here  were 
obtained  for  7^=463.15  K  and  Cbreak= 0.5  Cin. 
All  other  simulation  parameters  are  the  same  as 
in  the  simulations  of  adsorption  and  desorption. 
The  change  of  the  outlet  concentration  from  the 
adsorber  during  the  first  3  cycles  in  is  shown  in 
Figure  92,  the  maximal  temperature  in  the 
adsorbent  beds  in  Figure  93  and  the 
condensation  rate  in  Figure  94.  The  quasi¬ 
steady  state  is  practically  achieved  after  the 
second  cycle. 

It  should  be  noticed  that  the  TSA  cycle  in  the  2- 
cartridges  adsorber  lasts  considerably  longer 
that  in  the  1 -cartridge  adsorber.  This  is 
primarily  a  result  of  a  much  longer  adsorption 
half-cycle. 
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Figure  94.  Condensation  rate  in  the  2-cartridges 
adsorber  for  3  cycles 
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5.  Analysis  of  the  influence  of  process  parameters  and  optimization  of  the  TSA  process 

In  order  to  analyze  the  results  of  the  TSA  processes,  some  performance  criteria  should  be  defined,  as 
measures  of  the  process  quality.  We  define  the  following  criteria: 

•  Adsorption  time  (duration  of  the  adsorption  half-cycle),  xA 

•  Desorption  time  (duration  of  the  desorption  half-cycle),  xD 

•  Average  outlet  concentration  during  adsorption,  <Cout>  A 

•  Average  outlet  concentration  during  desorption,  <Cout>D 

•  Separation  factor  defined  as  the  ratio  of  the  average  outlet  concentrations  during  adsorption 
and  desorption,  as  =<c  „P>d/  <Cour>A 

•  Purification  factor,  defined  as  the  ratio  of  the  outlet  and  inlet  concentrations  during  adsorption, 

&p—  <"Goup>  a/  f^in,A 

•  The  total  amount  of  condensed  liquid  per  cycle,  Lcond 

•  Used  electric  energy  per  cycle,  Qd 

•  Used  energy  per  mol  of  condensate,  QeilLcond 

In  principle,  the  following  can  be  considered  as  favourable: 

Long  adsorption  time 

Short  desorption  times 

Low  outlet  concentration  during  adsorption 

High  outlet  concentrations  during  desorption 

High  separation  factor 

Low  purification  factor 

Large  amount  of  condensed  liquid  per  cycle 

Low  energy  spent  per  cycle 

Low  energy  spent  per  mole  of  condensate. 

In  this  section,  we  will  analyze  how  these  criteria  are  influenced  by  the  change  of  some  input 
parameters  of  the  TSA  process.  The  influence  of  the  following  parameters  was  chosen  for  analysis: 

The  maximal  adsorbent  temperature  at  which  the  TSA  cycle  is  switched  from  desorption 
to  adsorption  -  the  switch  temperature  Tsw; 

The  breakthrough  concentration  for  which  the  TSA  cycle  is  switched  from  adsorption  to 
desorption,  Cbreak; 

The  electrical  voltage  used  during  desorption  U0; 

The  inert  flow-rate  during  the  adsorption  part  of  the  cycle  GA; 

The  inert  flow-rate  during  the  desorption  part  of  the  cycle  GD. 

The  analysis  was  performed  by  changing  one  of  these  five  parameters  at  a  time  and  keeping  all  the 
others  constant.  The  “fixed”  parameters  have  the  same  values  as  in  the  simulations  presented  in  the 
previous  section.  For  each  combination  of  parameters  the  9  criteria  defined  previously  were 
calculated.  The  results  are  presented  in  tabular  form. 

The  results  presented  in  this  section  correspond  to  the  second  cycles  of  the  simulated  TSA  processes. 
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The  results  obtained  for  the  1 -cartridge  adsorber  are  presented  first. 


5.1.  Analysis  for  the  1-cartridge  adsorber 

In  Table  3  we  present  the  performance  criteria  chosen  for  analysis  obtained  for  three  different  switch 
temperatures  . 


Table  3.  Analysis  of  the  influence  of  the  switch  temperature  Tsw,  for  the  1 -cartridge  adsorber 


Criteria 

Switch  temperature  Tsw  [K]  I 

433.15 

463.15 

493.15 

<C0U^d  [mol/mol] 

0.00853 

0.00879 

0.00895 

<Cout>\  [mol/mol] 

2.65e-04 

2.00e-04 

1 ,84e-04 

as  D  /<Cout>A 

32.189 

43.950 

48.641 

®'P~<'C0ur> A  /Cjn  /\ 

0.265 

0.200 

0.184 

Lcond  [mol] 

0.0400 

0.0499 

0.0557 

Qel  [KJ] 

13.363 

16.735 

19.593 

Qel  !  LCOnd  [KJVlIlOl] 

334.07 

335.37 

351.76 

Td  [s] 

130 

158 

180 

TA  [si 

2370 

2440 

2600 

These  results  show  that  the  following  criteria  increase  with  increasing  the  switch  temperature:  the 
adsorption  and  desorption  times,  the  outlet  concentration  during  desorption  and  the  separation  factor, 
the  condensed  liquid  and  the  used  energy  per  cycle.  On  the  other  hand,  the  outlet  concentration  during 
adsorption  and  the  purification  factor  decrease  with  increasing  Tsw.  The  energy  consumption  per  mole 
of  condensate  is  nearly  constant.  It  could  be  concluded  that  the  increase  of  the  switch  temperature  is 
favorable  from  the  point  of  view  of  separation  and  purification  of  the  feed  stream,  as  well  as 
regeneration  of  the  adsorbed  vapor,  but  unfavorable  from  the  point  of  view  of  energy  consumption.  It 
should  also  be  remembered  that  the  switch  temperature  is  often  limited  by  the  temperature  of 
degradation  of  some  thermally  unstable  materials  used  in  the  process. 

In  Table  4  we  show  the  influence  of  the  breakthrough  concentration,  i.e.,  the  outlet  gas  concentration 
for  which  the  TSA  process  is  switched  from  adsorption  to  sorption. 


Table  4.  Analysis  of  the  influence  of  the  breakthrough  concentration  C break  f°r  the  1 -cartridge  adsorber 


Criteria 

Breakthrou; 

concentration  Chreuk  [mol/mol] 

2.5e-4  (25%  Cin) 

3.5e-4  (35%  Cin) 

5.0e-4  (50%  Cin) 

<C0Ut> D  [mol/mol] 

0.00835 

0.00854 

0.00879 

<C0ui> a  [mol/mol] 

1 ,28e-04 

1 .58e-04 

2.00e-04 

Ofs  ^ CffnP^X)  /<C0ut>A 

65.234 

54.051 

43.950 

A  /Cjn  /\ 

0.128 

0.158 

0.200 

Lcond  [mol] 

0.0397 

0.0437 

0.0499 

Qel  [KJ] 

15.268 

15.717 

16.735 

Qel  !  LCond  [KJVlIlOl] 

384.58 

359.66 

335.37 

Td  [s] 

143 

148 

158 

ta  [s] 

1800 

2050 

2440 

The  results  given  in  Table  4  show  that  the  increase  of  the  breakthrough  concentration  results  in  longer 
adsorption  and  desorption  times,  higher  outlet  concentrations  during  both  adsorption  and  desorption, 
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lower  separation  and  higher  purification  factors,  larger  amounts  of  condensed  liquid  and  energy  spent 
per  cycle,  but  lower  energy  consumption  per  mole  of  condensed  liquid. 


Table  5  shows  the  criterial  parameters  corresponding  to  simulations  obtained  with  3  different  values  of 
the  supply  voltage  during  desorption. 


Table  5.  Analysis  of  the  influence  of  the  voltage  used  for  desorption,  C/0,  for  the  1 -cartridge  adsorber 


Criteria 

Supply  voltage  t/0  [V]  | 

10 

11 

12 

<Cout> D  [mol/mol] 

0.00895 

0.00879 

0.00869 

<C„,„>A  [mol/mol] 

1 ,96e-04 

2.00e-04 

2.08e-04 

®'S~<'C0ur> D  /<C0ut>A 

45.663 

43.950 

41.779 

&p~<Cout> A  /Cin,A 

0.196 

0.200 

0.208 

Lcond  fmoll 

0.0500 

0.0499 

0.0510 

Qe,  [KJ] 

19.626 

16.735 

15.433 

Qei  /  Lcond  [KJ/mol] 

392.52 

335.37 

302.61 

Td  [si 

224 

158 

122 

Ta  [s] 

2460 

2440 

2430 

When  the  supply  voltage  during  desorption  increases,  slightly  shorter  adsorption  and  considerably 
shorter  desorption  times  are  obtained.  The  outlet  concentrations  during  adsorption  and  desorption  both 
increase  with  the  increase  of  £/0,  but  only  slightly,  as  well  as  the  purification  factor.  The  separation 
factor  decreases  somewhat,  and  the  used  energy  per  cycle  and  per  mole  of  condensate  considerably. 
The  amount  of  condensed  liquid  per  cycle  doesn’t  change  considerably. 

Tables  6  shows  the  influence  of  the  flow-rate  during  adsorption,  GA,  and  Table  7  of  the  flow-rate 
during  desorption,  GD,  on  the  criteria  chosen  for  comparison.  For  the  chosen  range  of  flow-rates  and 
for  the  parameters  used  for  simulation,  these  influences  are  not  significant. 


Table  6.  Analysis  of  the  influence  of  the  flow-rate  during  adsorption  GA,  for  the  1 -cartridge  adsorber 


Criteria 

Flow-rate  during  adsorption  GA  [mol/s] 

0.08 

0.09 

0.10 

<C0Ut> D  [mol/mol] 

0.00878 

0.0088 

0.00879 

<Cout>A  [mol/mol] 

2.02E-04 

2.01  E-04 

2.00e-04 

®'S~<'C0ut'> D  oup* \ 

44.7959 

43.7811 

43.950 

<'^our>\  /Cjn  /\ 

0.202 

0.201 

0.200 

Lcond  [mol] 

0.0502 

0.0506 

0.0499 

Qe,  [KJ] 

16.737 

16.842 

16.735 

Qei  /  Lcond  [KJ/mol] 

333.41 

332.85 

335.37 

Td  [s] 

158 

159 

158 

Ta  [s] 

2470 

2460 

2440 
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Table  7.  Analysis  of  the  influence  of  the  flow-rate  during  desorption  GD,  for  the  1 -cartridge  adsorber 


Criteria 

Flow-rate  during  desorption  GD  [mol/s]  [ 

0.01 

0.02 

0.03 

<C„ui>\)  [mol/mol] 

0.00879 

0.00868 

0.00886 

<C«r>A  [mol/mol] 

2.00e-04 

2.00e-04 

2.00e-04 

av  <Couf> D  A 

43.950 

44.8731 

43.9801 

Q'P~<'Cour> a  /Cjn  /\ 

0.200 

0.200 

0.200 

Lcond  [mol] 

0.0499 

0.0528 

0.0526 

Qel  [KJ] 

16.735 

17.471 

17.900 

Qel  !  LCOnd  [KJ/lHOl] 

335.37 

330.89 

340.30 

Td  [s] 

158 

165 

169 

XA  [si 

2440 

2410 

2430 

5.2.  Analysis  for  the  2-cartridges  adsorber 

As  our  analysis  for  the  1 -cartridge  adsorber  showed  that  the  flow-rates  don’t  influence  the  TSA 
process  much,  only  three  input  parameters  were  varied  for  the  2-cartridges  adsorber:  the  switch 
temperature,  the  breakthrough  concentration  and  the  electric  voltage.  Their  influence  on  the  analyzed 
performance  criteria  of  the  TSA  process  is  shown  in  Tables  8,  9  and  10,  respectively. 


Table  8.  Analysis  of  the  influence  of  the  switch  temperature  Tsw ,  for  the  2-cartridges  adsorber 


Criteria 

Switch  temperature  Tsw  [K]  [ 

433.15 

463.15 

485.28 

<C0Up> D  [mol/mol] 

0.00921 

0.00933 

0.00936 

<C0ui> a  [mol/mol] 

2.78E-04 

2.30E-04 

2.22E-04 

^s~<'^out> D  A 

33.1295 

40.4948 

42.1622 

&P  <(-'our> A  /Cin,A 

0.278 

0.230 

0.222 

Lcond  [mol] 

0.114 

0.1270 

0.131 

Qel  [KJ] 

34.648 

41.890 

44.665 

Qei  /  Lcond  [KJ/mol] 

303.93 

329.84 

340.95 

xD  [s] 

168 

198 

208 

Ta  [s] 

8090 

8570 

8720 

Table  9.  Analysis  of  the  influence  of  the  breakthrough  concentration  C break,  f°r  the  2-cartridges  adsorber 


Criteria 

Breakthrouj 

concentration  Cbreak  [mol/mol]  j 

2.5e-4  (25%  Cin) 

3.5e-4  (35%  Cin) 

5.0e-4  (50%  C,„) 

<C’(,«f>D  [mol/mol] 

0.00914 

0.00917 

0.00933 

<Cout>\  [mol/mol] 

1 .49E-04 

1 .92E-04 

2.30E-04 

as  <C()Ut>D  /<C()ut>A 

61 .3423 

47.7604 

40.4948 

®'P~<'C0ur> A  /Cjn  /\ 

0.149 

0.192 

0.230 

Lcond  [mol] 

0.0879 

0.1030 

0.1270 

Qel  [KJ] 

35.986 

37.531 

41.890 

Qei !  Lcond  [KJ/mol] 

409.40 

364.38 

329.84 

Xd  [s] 

170 

179 

198 

Xa  [si 

5720 

6840 

8570 
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Table  10.  Analysis  of  the  influence  of  the  voltage  supply  £/0,  for  the  2-cartridges  adsorber 


Criteria 

Voltage  supply  t/0  [V]  j 

10 

11 

12 

<C0Up> D  [mol/mol] 

0.00937 

0.00933 

0.00924 

<C0ui> a  [mol/mol] 

2.29E-04 

2.30E-04 

2.51  E-04 

^s~<'^out'> D  KCout>X 

40.9170 

40.4948 

36.8127 

&P  A  /Cjn  /\ 

0.229 

0.230 

0.251 

Lcond  [mol] 

0.1197 

0.1270 

0.1235 

Qe,  [KJ] 

45.837 

41.890 

35.899 

Qel  /  LCond  [KJ/lIlOl] 

382.93 

329.84 

290.68 

Td  [si 

265 

198 

143 

Ta  [s] 

8440 

8570 

8600 

Similar  conclusions  as  for  the  1 -cartridge  adsorber  can  be  withdrawn. 

Comparison  of  the  results  obtained  for  the  1 -cartridge  and  2-cartridges  adsorbers  shows  that  the 
adsorption  times  are  more  than  double  and  the  desorption  times  about  50%  longer  for  the  2-cartridges 
adsorber.  The  separation  and  the  purification  factors  are  similar  for  both  adsorbers.  On  the  other  hand, 
the  amount  of  condensed  liquid  and  the  used  energy  per  cycle  for  the  2-cartridges  adsorber  are  more 
than  double  than  for  the  1 -cartridge  adsorber.  Nevertheless,  the  values  of  the  energy  spent  per  mole  of 
condensate  are  similar  for  both  adsorber  types. 


5.3.  Discussion  of  the  results  and  some  recommendations  for  choosing  the  optimal  input 
parameters 

The  inspection  of  the  results  listed  in  Tables  3-10  shows  that  unexpectedly  low  values  of  the  amount 
of  condensed  liquid  per  cycle  are  obtained.  The  reason  for  that  could  be  in  the  fact  that  the  desorption 
step  in  our  model  is  ended  when  the  predefined  switch  temperature  is  reached,  while  the  condensation 
process  is  still  not  finished,  as  can  be  seen  from  Figure  53  (similar  picture  is  obtained  for  the  2- 
cartridges  adsorber).  After  that  moment,  the  simulation  switches  to  the  adsorption  model,  which 
doesn’t  assume  condensation.  In  order  to  resolve  this  imprecision,  an  additional  model  segment  should 
be  incorporated  in  the  integral  TSA  model,  for  an  additional  step  of  the  TSA  process  of  desorption 
with  in-vessel  condensation,  but  without  electrical  heating.  We  expect  that  such  a  model  would  predict 
larger,  more  realistic  amounts  of  condensate  per  cycle,  especially  for  the  cases  with  short  desorption 
times.  The  spent  energy  per  cycle  would  remain  unchanged,  so  we  expect  that  much  lower  values  of 
the  energy  consumption  per  mole  of  condensate  would  be  obtained. 

Regarding  the  influence  of  the  input  parameters  on  the  TSA  performance  criteria,  there  is  no  ideal 
solution,  i.e.  all  input  parameters  influence  some  of  the  criteria  in  a  favorable  and  some  in  an 
unfavorable  way.  For  that  reason,  the  recommendations  for  choosing  optimal  input  parameters  depend 
on  the  priorities  of  the  performance  criteria.  Here  are  some  examples: 

If  purification  of  the  gas  stream  is  the  main  priority,  lower  values  of  the  breakthrough 
concentration  should  be  used.  Higher  switch  temperature  and  lower  supply  voltage,  also 
result  with  more  efficient  purification,  to  some  extent. 

If  regeneration  of  the  vapor  in  the  liquid  form  is  the  priority  of  the  process,  one  should 
use  high  switch  temperature,  large  breakthrough  concentration  and  high  voltage. 

If  energy  efficiency  of  the  process  is  the  primary  goal,  low  switch  temperature,  large 
breakthrough  concentration  and  high  voltage  should  be  used. 
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6.  Conclusions  and  possibilities  for  model  improvement 


In  this  phase  of  the  project  we  have  developed  two  complex  models  for  simulation  of  integral  TSA 
cycles  in  adsorbers  with  one  or  two  annular,  radial-flow,  cartridge-type  adsorption  beds.  Each  model 
was  built  of  three  segments,  corresponding  to  three  steps  of  the  TSA  process:  adsorption, 
electrothermal  desorption  and  electrothermal  desorption  with  in-vessel  condensation.  The  individual 
segments  were  built  using  Femlab  software,  while  Femlab  with  Matlab  was  used  for  building  the 
integral  TSA  models. 

The  models  were  successfully  used  for  simulation  of  separate  stages  and  complete  TSA  cycles  in  the 
1 -cartridge  and  2-cartridges  adsorbers  and  for  analysis  of  the  influence  of  the  main  input  parameters 
on  the  system  performance.  Some  basic  recommendations  regarding  the  choice  of  the  input  parameters 
of  the  TSA  process  were  given  based  on  it,  in  order  to  show  that  the  developed  models  can  be  used  for 
optimization. 

The  models  were  developed  and  run  on  Pentium  4  PCs,  with  1  GB  RAM  memory.  During  their 
development,  convergence  problems  occurred  frequently.  Regarding  convergence,  the  models  are 
pretty  sensitive  on  the  choice  of  the  model  parameters.  The  computer  times  needed  for  simulations 
were  generally  long,  but  differed  from  case  to  case.  In  principle,  the  computer  times  needed  to  solve 
the  2-cartridges  model  were  much  longer  than  for  the  1 -cartridge  model. 


In  spite  of  the  great  complexity  of  the  developed  models,  there  are  still  some  possibilities  for  their 
upgrading  in  order  to  make  them  more  realistic.  Some  of  them  would  be: 

Adding  a  model  segment  defining  the  end  of  the  desorption  half-cycle  without  heating, 
between  switching  the  power  off  and  the  end  of  condensation. 

Adding  relations  for  calculation  of  some  physical  and  transport  parameters  and 
coefficients  as  functions  of  the  temperatures,  concentrations  and  velocities  in  the  system, 
instead  of  using  constant  values. 

Adding  a  model  for  the  heat  balance  of  the  adsorber  wall. 
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Nomenclature 


a  (m2/m3) 
b(  K'1) 

C  (mol/mol) 
C*(mol/mol) 

C break  (mol/lUOl) 

Csat  (mol/mol) 
cpg  (J/mol/K) 
cpl  (J/mol/K) 
cPs  (J/g/K) 
cpv  (J/mol/K) 

Dmr  (mol/cm/s) 
Dmz  (mol/cm/s) 
Ahg  (W/K/cm) 
Ahs  (W/K/cm) 

E  (J/mol) 
g  (cm/s2) 

G  (mol/s) 

A  (cm) 

/;,  (cm) 

/?2  (cm) 

Ab  (J/cm2/K) 
h8l  (J/cm2/K) 
hs2  (J/cm2/K) 
hWg  (J/cm2/K) 

J  (A/cm2) 

Jcond  (mol/cm2/s) 
k  (cm2) 

km  (mol/cm2/s) 

km  i  (mol/cm2/s) 

£m2  (mol/cm2/s) 
tube 


-  Specific  surface  area 

-  Temperature  coefficient  of  the  bed  electrical  resistivity 

-  Adsorbate  concentration  in  the  gas  phase 

-  Adsorbate  concentration  in  the  gas  phase  in  equilibrium  with  the  solid  phase 

-  Breakthrough  concentration 

-  Saturation  concentration 

-  Specific  heat  capacity  of  the  inert  gas 

-  Specific  heat  capacity  of  liquid  adsorbate 

-  Heat  capacity  of  the  solid  phase 

-  Specific  heat  capacity  of  gaseous  adsorbate 

-  Radial  mass  transfer  dispersion  coefficient 

-  Axial  mass  transfer  dispersion  coefficient 

-  Heat  diffusivity  of  the  gas  phase 

-  Heat  diffusivity  of  the  solid  phase 

-  Adsorption  energy  of  the  adsorbate  (D-R  eq.) 

-  Gravitation  constant 

-  Flow  rate  of  the  inert  gas 

-  Bed  axial  dimension  (Fig.  5) 

-  Adsorber  dimension  (Fig.  5) 

-  Adsorber  dimension  (Fig.  5) 

-  Solid  to  gas  heat  transfer  coefficient 

-  Heat  transfer  coefficient  from  the  solid  phase  to  the  gas  phase  in  the  central  tube(s) 

-  Heat  transfer  coefficient  from  the  solid  phase  to  the  gas  phase  in  the  annular  tube 

-  Gas  to  ambient  heat  transfer  coefficient  (heat  losses) 

-  Current  density 

-  Condensation  flux 

-  Asorbent  bed  permeability 

-  Mass  transfer  coefficient  in  the  adsorbent  bed 

-  Mass  transfer  coefficient  between  the  adsorbent  bed  and  the  gas  in  the  central  tube 

-  Mass  transfer  coefficient  between  the  adsorbent  bed  and  the  gas  in  the  annular 


Kond  (mol/s) 
Lcond  (mol) 

P  (Pa) 

Pa  (Pa) 


-  Flow-rate  of  the  condensed  liquid 

-  Total  amount  of  the  condensed  liquid 

-  Gas  pressure 

-  Adsorbate  partial  pressure 
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Pc  (Pa) 

-  Critical  pressure 

P  (Pa) 

-  Adsorbate  saturation  pressure 

Qel  (W) 

-  Rate  of  heat  generation  (equal  to  electric  power) 

Qe,(  J) 

-  Heat  generation  (equal  to  electric  energy) 

q  (mol/g) 

-  Adsorbate  concentration  in  the  solid  phase 

Rg  (J/mol/K) 

-  Gas  constant 

r  (cm) 

-  Radial  dimension 

r,  (cm) 

-  Adsorber  dimension  (Fig.  5) 

r2  (cm) 

-  Adsorber  dimension  (Fig.  5) 

r3  (cm) 

-  Adsorber  dimension  (Fig.  5) 

r4  (cm) 

-  Adsorber  dimension  (Fig.  5) 

Ta  (K) 

-  Ambient  temperature 

T’g(K) 

-  Gas  phase  temperature 

Tc(  K) 

-  Critical  temperature 

7r(K) 

-  Referent  temperature 

T’s  (K) 

-  Solid  phase  temperature 

Ts,  (K) 

-  Switch  temperature  (from  desorption  to  adsorption) 

Tw  (K) 

-  Wall  temperature 

*(s) 

-  Time 

C/(V) 

-  Electric  potential 

t/o(V) 

-  Supply  voltage 

w  (cm/s) 

-  Radial  superficial  gas  velocity 

EP,,  FP* 

VPD  -  Wagner  constants 

v  (cm/s) 

-  Axial  superficial  gas  velocity 

ff(cm3/g) 

-  Adsorbate  volume  per  1  g  of  the  adsorbent  mass 

(cm3/g) 

-  Total  volume  of  micropores  (D-R  equation) 

Greek  letters: 

ap 

-  Purification  factor 

ar 

-  Regeneration  factor 

as 

-  Separation  factor 

A/Ads  (J/mol) 

-  Heat  of  adsorption 

Aflcond  (J/mol) 

-  Heat  of  condensation 

£b 

-  Bed  porosity 

p  (Pa/s) 

-  Dynamic  viscosity 

p  (Qcm) 

-  Electric  resistivity 

po  (Qcm) 

-  Electric  resistivity  at  referent  temperature  TR 
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pg  (mol/cm3) 

-  Gas  phase  density 

pb  (g/cm3) 

-  Adsorbent  bed  density 

Pa  (g/cm3) 

-  Adsorbate  density 

x(s) 

-  “Half’  cycle 

Subscripts: 

A 

-  adsorption 

b 

-  bed 

D 

-  desorption 

g 

-  gas 

in 

-  inlet 

out 

-  outlet 

P 

-  previous  (initial) 

r 

-  radial 

s 

-  solid  phase 

it 

-  inner  (central)  tube 

ot 

-  outer  (annular)  tube 

z 

-  axial 

Special  characters: 

<>  -  time  average  value 
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